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ABSTRACT 


In 1979, W. B. Colson and S. K. Ride proposed a new kind of electron accelerator using a 
uniform magnetic field in combination with a circularly-polarized laser field. A key con¬ 
cept is to couple the oscillating electric field to the electron’s motion so that acceleration 
is sustained. This dissertation investigates the performance of the proposed laser accelera¬ 
tor using modern high powered lasers and magnetic fields that are significantly improved 
over the values in the original analysis. The relativistic equations of motion for the elec¬ 
tron are solved analytically and for the first time numerically for more general results than 
previously reported. Realistic beam effects are investigated to determine the constraints 
on initial electron energy, injection angle, and phase to achieve optimal performance. This 
analysis indicates that the selection of an electron’s initial conditions are not as stringent 
as previously thought. Simulations indicate a university scale system can achieve electron 
energy gain of over one GeV in a meter. 


v 



THIS PAGE INTENTIONALLY LEFT BLANK 


vi 



Table of Contents 


1 Background 1 

1.1 Motivation and Research Questions. 1 

1.2 History of Particle Accelerators. 2 

2 Laser Accelerator 9 

2.1 Laser Accelerator System Overview. 9 

2.2 Laser Source. 9 

2.3 Electron Injector. 14 

2.4 External Magnetic Field. 18 

3 Equations of Motion 23 

3.1 Relativistic Lorentz Force. 23 

3.2 Equations of Motion. 25 

3.3 Integral Representation. 31 

3.4 Curtate Cycloids. 33 

3.5 Resonance Condition. 37 

3.6 Electron Phase. 38 

3.7 Maximum Energy Gain. 39 

4 Numerical Methodology 41 

4.1 Finite-difference Method. 41 

4.2 Initial Conditions on Resonance (Aft) = 0). 42 

4.3 Electron Trajectory. 44 

4.4 Synchrotron Radiation. 45 

4.5 Numerical Integration. 47 

5 Numerical Analysis Results 49 

5.1 On-resonance with 0 o — n . 49 

5.2 Electron Phase Performance. 53 

vii 
























5.3 Off-resonance Performance. 60 

5.4 Magnetic Field Effects. 66 

5.5 Diffraction Effects. 67 

6 Conclusions 73 

Appendix: Simulation Codes 75 

A.l Laser Accelerator Simulation Core. 75 

A.2 Laser Accelerator Matlab Graphics Script. 92 

List of References 99 

Initial Distribution List 103 


viii 








List of Figures 


Figure 1.1 Cockroft-Walton Cascade Generator . 3 

Figure 1.2 Van de Graaff Accelerator. 4 

Figure 1.3 Linear Accelerator. 5 

Figure 1.4 Cyclotron Accelerator. 6 

Figure 1.5 Livingston Chart. 7 

Figure 2.1 Chirped Gaussian Pulse. 11 

Figure 2.2 Diagram Illustrating Chirped Pulse Amplification . 12 

Figure 2.3 Texas Petawatt Laser System Layout. 14 

Figure 2.4 Cross-section of a Triple Spoke Cavity. 16 

Figure 2.5 Naval Postgraduate School SCRF Gun. 17 

Figure 2.6 Cutaway Diagram of Hybrid Magnet at the National High Magnetic 

Field Laboratory. 19 

Figure 2.7 Cutaway Diagram of a Pulsed Magnet System at the National High 

Magnetic Field Laboratory . 21 

Figure 3.1 Schematic of Laser Accelerator. 27 

Figure 3.2 Standard, Curtate, and Prolate Cyloids . 34 

Figure 3.3 Diagram of Initial Electron Distribution in the Longitudinal Direc¬ 
tion . 38 

Figure 4.1 Injection Angle Required for Resonance. 43 

Figure 4.2 Allowed Values of Injection Angles. 44 

Figure 5.1 Comparison of Analytic Solution to Numerical Simulation Results 

of an Electron in the Laser Accelerator. 50 






















Figure 5.2 Numerical Simulation Results For Two Different Electrons in the 

Laser Accelerator for 45 T and 100 T External Magnetic Fields . 51 

Figure 5.3 Electron Trajectories in the Laser Accelerator for Two Different Ex¬ 
ternal Magnetic Field Cases. 53 

Figure 5.4 Numerical Simulation of Laser Accelerator for Various Initial 9 0 and 

y 0 = 130 . 55 

Figure 5.5 Numerical Simulation of Laser Accelerator for Various Initial 9 0 and 

y 0 = 125 . 56 

Figure 5.6 Numerical Simulation of Laser Accelerator for Various Initial 0 o and 

Yo =120 . 57 

Figure 5.7 Numerical Simulation of Laser Accelerator for Various Initial 0 o and 

Yo =119 . 58 

Figure 5.8 Numerical Simulation of Laser Accelerator for Minimum Energy 

Injection for Resonance. 59 

Figure 5.9 Numerical Simulation of Laser Accelerator Over All 0 o and Small 

Range of Injection Angles. 61 

Figure 5.10 Numerical Simulation of Laser Accelerator Over All 9 a and Small 

Range of Injection Angles. 62 

Figure 5.11 Numerical Simulation of Laser Accelerator Over All 9 a and Small 

Range of Injection Angles. 63 

Figure 5.12 Numerical Simulation of Laser Accelerator Over All 9 a and Small 

Range of Injection Angles. 64 

Figure 5.13 Numerical Simulation of Laser Accelerator Over B„ and Initial En¬ 
ergy Yo . 65 

Figure 5.14 Magnetic Field Variation Along z with Maximum Field at z w = 0.5 

m. 67 

Figure 5.15 Magnetic Field Variation Diagram Defines 9\ and 02 for Calculating 

the Magnetic Field Along the Longitudinal Axis of a Solenoid . . 67 

Figure 5.16 Numerical Simulation Results of a Single Electron in the Laser Ac¬ 
celerator System with Magnetic Field Variation. 68 


x 
















Figure 5.17 


Laser Phase Evolution 


70 


xi 




THIS PAGE INTENTIONALLY LEFT BLANK 



List of Tables 

Table 2.1 Petawatt Laser System Specifications, after [2] 10 

Table 2.2 Available Magnets at the National High Magnetic Field Laboratory, 

after [18] 18 

Table 2.3 Available Pulsed Magnets at PFF LANL. SP = Short Pulse. MP = 

Mid-Pulse, after [19]. 20 

Table 2.4 Magnet Specification Summary, after [19]. 20 

Table 3.1 Equations and Constants of the Motion. 31 

Table 5.1 Input Parameters for Resonance Laser Accelerator Simulation . . 50 

Table 5.2 Results for Single Electron Simulation of Laser Accelerator at 45 T 

and 100 T . . .. 53 


xiii 










THIS PAGE INTENTIONALLY LEFT BLANK 


xiv 



List of Acronyms and Abbreviations 


a injection angle 
A ampere 

a third constant of the motion 

b fourth constant of the motion 

B„ external magnetic field strength 

/3 electron velocity relative to the speed of light 

/3.v x component of electron velocity relative to the speed of light 

j8y y component of electron velocity relative to the speed of light 

/3- z component of electron velocity relative to the speed of light 

/3j_ perpendicular component of electron velocity relative to the speed of light 

c speed of light 

CPA chirped pulse amplification 
C coulomb 
DC direct current 
e normalized emittance 
e first constant of the motion 
E, laser electric field amplitude 
k e electron wavenumber 
k r laser radiation wavenumber 
E c electric field amplitude 
eV electron volt 
fs femtosecond 
a p Gaussian pulse length 
Hz Hertz 

(Qj instantaneous angular frequency 

J joule 

K Kelvin 

kA kiloampere 

keV kilo-electronvolt 

A,- laser wavelength 


xv 




bp linear chirp parameter 
7 Lorentz factor 
m meter 

NHMFL National High Magnetic Field Laboratory 
rj second constant of the motion 
Nd neodymium 

NCRF normal conducting radio-frequency 
OPCPA optical parametric chirped pulse amplification 
//' momentum four-vector 
(f> azimuthal angle between x-axis and /3_ 

PW petawatt 

R radius of curvature of wavefronts 

0 angle between perpendicular velocity and laser electric field vector 
TPL Texas Petawatt Laser 
w optical mode waist 
A ft) resonance parameter 
( o c relativistic cyclotron frequency 

(Oe electric field analogue to relativistic cyclotron frequency 

C0 o central Gaussian wavelength 

(0 p pump laser angular frequency 

(O s signal laser angular frequency 

(O rs residual signal angular frequency 

ft),, laser radiation angular frequency 

£ electron phase 

Z, Rayleigh range 


xvr 



Acknowledgments 


I would like to express my deepest and utmost gratitude to my advisor, Prof. Bill Colson. 
I am grateful for his mentorship, patience, and friendship over the past ten years. I look 
forward to continued collaboration and discussion in my next tour of duty and beyond. 

I would like to thank Prof. Cohn and Prof. Blau for their guidance especially in helping 
me reduce my long-winded sentences to something clear and concise. I would l ik e to thank 
Prof. Newman, Prof. Karunasiri, and Prof. Luscombe for their support. I will miss Prof. 
Luscombe’s classroom discussions the Grateful Dead, Monty Python, and random movie 
quotes. 

I would also like to thank Prof. Giraldo, Prof. Frenzen, and my cousin Prof. Dan Gefroh 
for their valuable discussions about this dissertation. 




THIS PAGE INTENTIONALLY LEFT BLANK 


xviii 



CHAPTER 1: 
Background 


In 1979, W. B. Colson and S. K. Ride [1] proposed a new kind of accelerator using a uni¬ 
form, magnetic field in combination with a circularly-polarized laser field. Lasers can have 
strong electric fields with frequencies high enough to avoid corona formation and break¬ 
down. The key is to couple the electric field to the electron such that it always “sees” a 
relatively steady field so that acceleration is sustained. Colson and Ride [1] suggested that 
if the external magnetic field is aligned along the light-propagation axis and the magnetic 
field strength, laser frequency, and initial electron velocity are selected judiciously, the 
composite field forces are resonant indefinitely, so that the electron continues to be acceler¬ 
ated over the interaction length. For laboratory magnetic fields of 45 T and a petawatt laser, 
it will be shown through numerical simulations, electrons could be accelerated from about 
100 million-electron-volts to about 1 giga-electron-volts (GeV) in about a meter. This is 
significantly better than considered in the original research paper [1], 

1.1 Motivation and Research Questions 

In today’s high powered lasers, electric fields on the order of 10 11 V/m are possible [2], 
Electric fields are transverse to propagation and oscillating at high frequencies on the order 
of 10 15 seconds [2], [3], This high frequency conspires to make the coupling of the elec¬ 
trons to the light difficult. Electromagnetic forces would average to zero due to the rapidly 
oscillating light field. Using an external uniform, static magnetic field on the order of 45 
T, would help provide sustained coupling between the electron and the alternating fields of 
the laser. This motivates a several questions that this dissertation will attempt to answer. 

1. Can modern lasers and magnetic fields provide increased coupling compared to Col¬ 
son’s and Ride’s original results? 

2. What are the engineering considerations for this laser accelerator system? 

3. How is success measured? 

4. What are the constraints, if any, for optimal performance? 
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5. How would real beam effects, such as emittance, affect performance? 


1.2 History of Particle Accelerators 

Particle accelerators have been used since the early 1900s to study the structure of matter. 
Since then, particle accelerators have also been used in applied research for industrial, med¬ 
ical, and military applications. Particle accelerators are broadly classified into two types: 
electrostatic or direct current (DC) accelerators and oscillating field or radio-frequency 
(RF) accelerators. 

An important discovery in the development of particle accelerators can be traced back to the 
work of E. Rutherford [4]. In 1906, Rutherford developed the theory of atomic scattering 
through his experiments with a-particles and atomic disintegration [4]. His subsequent 
gold-foil experiments with H. Geiger and E. Marsden put to rest the “plum-pudding” model 
of the atom in favor of the nucleus model. The Geiger-Marsden experiments proved that 
unaccelerated a-particles, produced by the spontaneous decay of radium atoms, could be 
scattered by over 90 degrees when fired at a gold foil. Scattering angles of this magnitude 
could only occur if a dense nucleus existed at each gold atom site [4]. 

By 1919, Rutherford believed he required on the order of a few million electron-volt (MeV) 
a-particles to split the lithium nucleus. This value exceeded the electrostatic accelerator 
technology at that time so research in this effort stalled [5]. In 1928, G. Gamov predicted, 
through quantum tunneling effects, that approximately 450 keV a-particles would suf¬ 
fice in splitting the lithium nucleus. Encouraged by Rutherford, Cockroft and Walton de¬ 
veloped a DC accelerator, shown in Fig. 1.1, to achieve this energy requirement. The 
Cockroft-Walton accelerator is a type of accelerator whereby static electric fields are used 
to accelerate charged particles [6]. Their design uses a voltage multiplier comprised of 
capacitors and diodes to achieve high DC voltage from a low DC or AC voltage source. 
By 1932, their accelerator reached 400 keV and was able to split the lithium nucleus. They 
were awarded the Nobel Prize in Physics in 1951 for this work. In DC or electrostatic 
accelerators, voltages of several MV can be reached today [5] . 

Another type of electrostatic accelerator is the Van de Graaff generator developed by R. J. 
Van De Graaff in the early 1930s, see Fig. 1.2. In this design, a moving belt is used to 
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transfer charge to a hollow metal sphere. Voltages of over 10 MV can be reached if used 
with an insulating gas. Corona discharge limits all electrostatic accelerators to voltages 
of several tens of MV due to the formation of ions in the device. The stripped electrons 
then collide and cause an avalanche of more electrons leading to a spark discharge in the 
accelerator and thus a breakdown of the electrostatic field [6], [7]. 



Figure 1.1: Cockroft-Walton cascade generator used to split lithium atom, from [7]. 
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Figure 1.2: Schematic of a Van de Graaff accelerator, from [7]. 


To overcome corona formation, G. Ising proposed in 1924 the use of radio-frequency (RF) 
fields [4]. Corona formation is avoided since the alternating voltage seen by molecules is of 
a high enough frequency that the electromagnetic force time averages to zero. In 1928, R. 
Wilderoe demonstrated this principle with 50 keV potassium ions [4]. Wilderoe’s device is 
a linear RF accelerator using drift tubes, as shown in Fig. 1.3. The RF source is connected 
to each drift tube so that an electron only feels the accelerating force of the electric field 
at each drift tube gap. Then, while the electron is traveling at constant velocity in the drift 
tube, the fields oscillate such that the electron is again accelerated in the next gap. Drift 
tubes have been replaced by metal cavities since the mid-1930s and by superconducting 
RF cavities since the 1960s [4]. Superconducting RF cavities minimize the heat losses in 
the cavity relative to the energy stored in the resonator achieving a very high Q-factor, and 
therefore high fields [4]. 
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Figure 1.3: Schematic of a linear accelerator, from [7]. 


The cyclotron, developed in the 1930s by E. O. Lawrence, makes use of a similar accel¬ 
erating structure, but is more compact due its design as shown in Fig. 1.4. Particles orbit 
in a plane between two poles of an external magnet at a frequency known as the cyclotron 
frequency. The external magnet provides a static magnetic field that holds the particles in 
the plane while the varying RF fields on each of the D-shaped electrodes, called “dees,” 
provide the accelerating force. The charged particles accelerate outward from the center of 
device in a spiral fashion. Due to the electron’s small mass, the cyclotron was limited in 
energy by relativistic effects [4] and was better suited for heavier particles or nuclei. 

A special type of cyclotron, called the “synchrocyclotron,” was developed by E. McMillian 
in the 1940s. In this design, the RF fields are varied to compensate for the relativistic 
effects. Wilderoe suggested a new design such that charged particles form the secondary 
loop of a transformer [4]. In this new device, called a “betatron,” the magnetic field is 
increased as the particles accelerate. The particle trajectories remain circular and with a 
constant radius [8]. Independently, D. Kerst built the first betatron in 1940 that produced 
2.2 MeV electrons. By 1950, betatron achieved energies of 300 MeV. 

To achieve higher energies and to maintain a constant radius, the magnetic field had to be 
synchronously adjusted with energy of the particle [7]. While a classical cyclotron uses 
both a constant guiding magnetic field and a constant-frequency electromagnetic field, its 
successor, the “isochronous cyclotron,” works by local variations of the guiding magnetic 
field, compensating for the increasing relativistic effects of the particles during accelera¬ 
tion. In a synchrotron, this adaptation is done by variation of the magnetic field strength 
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in time, rather than in space [6]. In proton synchrotrons, energies of over 1 TeV can be 
reached. Some energy is lost through synchrotron radiation that is emitted when charged 
particles are accelerated radially. 

During the latter half of the last century, particle accelerator energies nearly doubled every 
seven years. Since most exotic particles of interest do not exist freely, through consequence 
of Einstein’s famous equation relating energy and matter, exotic particles must be created 
through collisions at high-energies. The quest for higher energies to unlock new exotic 
particles has pushed accelerator technology over the last century. The exponential growth 
in accelerator energy can be summarized in a Livingston chart (see Fig. 1.5). Since the 
1980s, the pace of accelerator energies has fallen off. Today, the Large Hadron Collider is 
the most powerful accelerator in the world. Powered by a proton synchrotron beam line, it 
can reach energies of seven TeV. 

New types of devices using plasma waves to accelerate electrons are being developed. 
These include the laser wakefield accelerator, plasma wakefield accelerator, plasma beat- 
wave accelerator, and self-modulated laser wakefield accelerator. These technologies may 
have the potential to accelerate electrons to a few GeV in a about a meter [9]. 



Figure 1.4: Schematic of a cyclotron accelerator, from [7], 
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Figure 1.5: Exponential growth of accelerator beam energy, from [5] 
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CHAPTER 2: 
Laser Accelerator 


This chapter will discuss the proposed components for a laser accelerator system using real 
world examples to frame an engineering estimate for power, size, and cost estimates. 

2.1 Laser Accelerator System Overview 

This laser accelerator system has three major components: 

1. Laser source - provides the high electric field used in accelerating the electrons 

2. Electron injector - produces a beam of electrons and directs into the interaction region 
at an initial angle and energy 

3. External magnet - provides the magnetic field component of the Lorentz force that 
enables coupling between the electrons and laser light 

2.2 Laser Source 

High peak power lasers with large electric fields can be laboratory scale and provide the 
force for significant particle acceleration. As an example of high peak power laser tech¬ 
nology, the Texas Petawatt Laser (TPL) at the University of Texas at Austin’s Center for 
High Energy Density Science can reach energies of 200 J in 170 femto-seconds, or over 
1.1 petawatts (PW) peak power [2], [10], [11] . The laser system provides two beam lines 
directed into three different target chambers for scientific experiments. The main beam line 
is based on optical parametric chirped pulse amplification (OPCPA) followed by power am¬ 
plification in two types of neodymium-doped glass [2]. Several experiments planned for the 
Texas Petawatt Laser (TPL) include: production of fusion neutrons, study of laser heated 
matter, dynamics of shocked materials, radiative hydrodynamics, bright x-ray research, and 
particle acceleration [10]. 

TPL produces 170 fs pulses at 1057 nm central wavelength [11] . Users can select several 
energy levels from 60-200 J per pulse at a rate of one shot per hour, or 7 J per pulse at four 
shots per hour [11] . Table 2.1 summarizes the laser system specifications [2]. 
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Table 2.1: Petawatt Laser System Specifications, after [2] 



Short Pulse 

Long Pulse 

Pulse Energy 

Pulse Duration 

Wavelength 
Repetition rate 
Focused intensity 

200 J 

170 fs 

1057 nm 

1 shot/hour 

> 10 21 W/cm 2 

500 J 

2-20 ns 

527 nm 

1 shot/hour 

10 17 W/cm 2 


At the heart of the petawatt laser system are three OPCPA stages that are used to amplify 
the seed oscillator energy from mJ to nearly 200 J per shot. OPCPA techniques are the latest 
generation of chirp pulse amplification (CPA) technology. CPA has enabled scientists to 
push the envelope of peak power systems from gigawatts to petawatts in the past 30 years. 
CPA was first demonstrated by D. Strickland and G. Mourou in 1985 [3]. We will briefly 
describe how CPA and OPCPA work. 

Consider a Gaussian pulse written in the following form [3]: 

E{t) = £>-vV(<*V+V 2 ) (2.1) 


where E„ is the peak electric field strength, a p is a measure of the Gaussian pulse length, b p 
is the linear chirp parameter, and co () is the light’s central angular frequency. The first ex¬ 
ponential factor describes the Gaussian envelope and the second term describes the phase. 
The total instantaneous phase is given by [3]: 


0total — T bpt ■ 


The instantaneous frequency is defined as [3] 


(Oi{t) 


d ^total 

dt 


Therefore, the instantaneous frequency for the pulse given by Eq. (2.1) is [3] 


( 2 . 2 ) 


(2.3) 


O)i(t) = (0 o + 2b p t. 


(2.4) 
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As can be seen from Eq. (2.4), the instantaneous signal frequency increases with time that 
is called an up-chirp for b p > 0. Down-chirp pulses have their signal frequency decreasing 
with time with b p < 0. Exponential chirps are also possible. The pulse can acquire a chirp 
during transmission through a dispersive or nonlinear medium. Fig. 2.1 shows a chirped 
Gaussian pulse. 


E/E 0 


time 


Figure 2.1: Diagram illustrating chirped Gaussian laser pulse. The blue curve represents the highly 
chirped signal. The red curve represents the Gaussian envelope. 



Fig. 2.2 illustrates basic the concept of chirped pulse amplification [12]. An oscillator 
produces a short pulse (1), which is then fed into gratings that stretch the pulse temporally 
out by 10 3 - 10 5 times (2). The gratings cause the low-frequency components to travel a 
shorter path than the high-frequency components, thereby giving it an up-chirp. The result 
is a temporally stretched pulse with a low enough peak power to safely pass through the 
main laser gain medium usually Nd:glass or Tfisapphire (3). The amplified pulse is then 
recompressed with gratings giving it a down-chirp resulting in a very short and high energy 
pulse (4). 
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1) Oscillator produces short pulses 



2) Using gratings, the pulse 
Is stretched by 10,000 x 


3) Long pulses 
are now sate to 
amp Illy In laser 
amplifiers 


4) After amplification, 
the pulses are 
recompressed vtrith 
gratings 


5) Resulting pulse Is 
short <20-1000 fs) 
and high energy (> 1 J) 


Figure 2.2: Diagram illustrating chirped pulse amplification, from [12]. 


Parametric amplification is an effect where the dielectric polarization responds nonlinearly 
to the electric field of input light. A higher frequency laser, known as the pump, is shined 
onto the parametric crystal resulting in a splitting of the pump photon into a pair of lower 
frequency photons. The input consists of a pump and signal photons while the output 
consists of residual pump, signal, and idler photons respectively such that ft),- = (O rp — (O s . 
The phase relationship between signal and pump determines the direction of energy flow 
whether there is amplification or deamplification of the signal [13]. Several benefits of para¬ 
metric amplification include: a much wider range of acceptance wavelengths compared to 
laser amplifiers, higher gain per unit length of gain material in pulsed mode, good thermal 
properties since there is very little absorption in the nonlinear crystal. Disadvantages of 
parametric amplification include no pump energy acculumation in the gain medium, pre¬ 
cise pump and signal synchronization, and losses introduced by idler photons [14]. 

OPCPA combines the techniques of CPA and parametric amplification to achieve very high 
peak power in ultrashort pulses. The main thrust of OPCPA is to amplify a short laser 
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pulse while mitigating some deleterious nonlinear effects, such as self-focusing and phase 
modulation. TPL makes extensive use of OPCPA technology and begins with a 100 fs 
tunable oscillator [2]. The pulse is stretched to 2 ns as described by the CPA techniques 
above. There are three stages of parametric amplification using nonlinear crystals. Stage 
one consists of three beta barium-borate (BBO) crystals and a Nd:YAG pump laser. Stage 
two uses a pair of BBO crystals also pumped by a Nd:YAG laser. Stage three uses a 532 
nm, 4 J, custom laser pump and two potassium deuterium phosphate (KDP) crystals. The 
final preamplifier stage is a Ti:sapphire crystal. The light source for pumping is residual 
light from the YAG in the 2nd or 3rd stage. This then feeds the main amplifier comprised 
of two Nd:glass amplifiers. Pulse energy is increased to about 300 J. Finally, light is passed 
to the compressor where it is housed in a 7x11 foot vacuum vessel. A multilayer dielectric 
diffraction grating brings the pulse energy down to approximately 200 J, 150 fs and is 
directed to the target bay for experiments [2]. 

OPCPA enables high gain, small thermal management issues (due to pulse stretching), 
flexible peak gain, and spectral gain [10]. 

TPL is electrically powered with a large capacitor bank next to the laser/target bay rooms. 
Total stored energy is approximately 2 MJ. Fig. 2.3 illustrates the layout and size of system 
with the laser bay and target area comprising of approximately 1500 square meters [10]. 
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Figure 2.3: A schematic diagram that illustrates laser bay, target bay, and capacitor bank rooms, 
from [10] 


2.3 Electron Injector 

Electron injectors, or electron guns, come in many different varieties and configurations. 
The basic components of all electron injectors are the cathode for emitting electrons into 
free space, a beam focusing structure usually comprised of electromagnets, and the subse¬ 
quent accelerating cavities. Injectors can be classified by their: 

1. Emission mechanism: thermionic, photocathode, or field emission cathode. 

2. Electric field generation: static field or direct current (DC) guns, radio-frequency 
(RF) guns, normal conducting RF (NCRF) guns, or superconducting RF (SCRF) 
guns. 

Other performance properties that are important include: average beam current, peak beam 
current, electron pulse length, emittance, physical size, cooling, reliability, and cost. 
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Thermionic emission of electrons uses the heat-induced ejection of electrons from a metal 
or non-metal surfaces. A common example of this is a vacuum tube where a hot filament is 
used to impart sufficient kinetic energy to electrons so they can overcome the metal’s work 
function at the surface. A positive anode is placed some distance away from the surface 
so that the resulting electrostatic force will pull the electrons into vacuum away from the 
cathode where they can be further accelerated. 

Photocathode emitters eject electrons from the cathode by the photoelectric effect. A drive 
laser source shines light energy onto a metal or semiconductor surface [15]. The quantum 
efficiency (QE) of a cathode is the ratio of the number of electrons generated to the number 
of incident photons. Metals, such as copper, tend to have a low QE and require a UV drive 
laser but have a fast response time. Semiconductors, such as gallium arsenide or cesium 
telluride, have higher QE, and only require a visible drive laser, but have a slower response 
time. The UV drive laser results in a fast response time while the visible drive laser gives a 
slower response time. 

Field emission cathodes eject electrons from a surface due to a nearby high electric field 
gradient that lowers the work function so that electrons can tunnel through the barrier and 
escape. Field emission cathodes are also considered cold cathodes since a heating element 
is not present. Field emitters structures tend to have spikes or sharp tips that greatly enhance 
the electric field at the tip near the surface. Recent papers have shown that carbon nanotubes 
may be excellent field emitters [16]. 

Following the cathode emission section of the electron gun, DC or RF fields provide the 
initial stage of acceleration. RF cavity shapes can vary from pillbox, elliptical, quarter- 
wave, or spoked. Spoke cavities, see Fig. 2.4, are a coaxial half-wave length cavity with an 
outer conductor turned ninety degrees so that its axes is directed along the beam path [17]. 

Electron guns with a photocathode can be combined with a DC, NCRF, or SCRF acceler¬ 
ating cavity section. NCRF gun operations span a broad range of RF frequencies ranging 
from a few MHz to several GHz. The RF field and the drive laser pulses must be care¬ 
fully synchronized so that the electrons ejected via the photoelectric effect are properly 
accelerated. 
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Figure 2.4: Cross-section of a triple spoke cavity, from [17] 


In NCRF, resistive losses cause the cavity walls to heat up and require cooling systems 
based on a water heat-exchanger. Heating is greatly reduced through use of a supercon¬ 
ducting cavity but system complexity and cost is increased. The liquid-He requirement for 
some super-conducting cavities may be somewhat mitigated by choosing a lower RF band. 
For RF of less than about 500 MHz, a liquid-He system of only 4K is sufficient. A 2K 
liquid-He system is required for RF greater than about 500 MHz. 

Electron guns can be continuous wave or pulsed. A particle synchronized with the RF field 
is called a synchronous particle. Particles bunched near the synchronous particle in the 
oscillation of the RF field is called a bucket. Electrons can be present in each RF bucket, 
or in pulsed guns where certain RF buckets are filled. Typical frequencies are a few MHz 
to several gigahertz with bucket charges from pC to nC per bunch. 

As mentioned earlier, electron beam emittance is an important quality to consider. It is 
a measure of particle spread in position and angle and is typically measured in units of 
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mm-mrad. Normalized emittance, £ n is defined as 


£ n = Yo<r><e> (2.5) 

where y 0 is the beam Lorentz factor, < r > is the average spread in particle position and 
<0> is the average spread in particle angle. As we will see later, the electron beam positions 
must be well inside the optical mode of the laser to be accelerated. For a relatively good 
beam emittance of 10 mm-mrad and y () — 130, the spread in position and angles must be 
less than 0.08 mm-mrad. We shall see later that this is a sufficiently small spread in angles 
in order to achieve optimal energy gain. 

The Naval Postgraduate School and Niowave, Inc have designed and developed the Mark I 
SCRF injector (see Fig. 2.5). This design uses a quarter wave structure, which are known 
to achieve a high accelerating efficiency. Furthermore, operating at a low RF band of 
350 MHz allows use of a 4 K helium refrigeration system thereby greatly minimizing the 
complexity of total system. The Mark I uses an RF-coupler with a nominal electron beam 
energy of a few MeV. 



Figure 2.5: Naval Postgraduate School SCRF Gun 
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2.4 External Magnetic Field 

The external magnetic field in the laser accelerator system serves to provide coupling be¬ 
tween the electron beam trajectory and the laser’s transverse electric field. It will be shown 
that magnetic field strengths on the order of 30 T to 60 T provide for sufficient electron 
coupling. We will find that higher magnetic fields allow the initial electron energy and 
electron injection angle to be much lower. 

An example of high magnetic field technology, consider the National High Magnetic Field 
Laboratory (NHMFL). NHMFL currently holds the record for the strongest continuous 
static field in the world at 45 T. NHMFL is a federal-state research facility in partnership 
with Florida State University, the University of Florida, and Los Alamos National Labo¬ 
ratory under a National Science Foundation grant. Its mission: “To provide the highest 
magnetic fields and necessary services for scientific research conducted by users from a 
wide range of disciplines, including physics, chemistry, materials science, engineering, bi¬ 
ology and geology” [20]. There are several magnet types available for scientists to conduct 
basic research. Table 2.2 summarizes some of NHMFL’s magnets that are available [18]. 

Table 2.2: Available Magnets at the National High Magnetic Field Laboratory, after [18] 


Field [T] 

Bore [mm] 

Power [MW] 

Resistive and Hybrid Magnets 



20 

195 

20 

24.5 

32 

14 

25 

52 

19 

27 

32 

14 

30 

32 

16 

33 

32 

16 

45 

32 

15 

Superconducting Magnets 


Temperature 

15 

45 

lOmK - 1 K 

17.5/19.5 

52 

0.4 - 300 K 

18/20 

52 

20 mK - 2 K 


The Hybrid Magnet at the NHMFL is the world record holder for continuous magnetic 
field (see Fig. 2.6). It consists of a room temperature magnet (31 T) surrounded by a 
superconducting magnet (14 T) to produce a continuous, static, 45 T field in its 32 mm 
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bore over one m length. The power supply is rated for 40 MW [20]. The cooling system 
is buffered by a 3800 m 3 cold water storage tank that is maintained at 7°C by a 28 MW 
chiller [20], [21]. 



Electrical and 
Helium Supply tor 
Superconducting 
Magnet 


Outsert Magnet 
Cryostat 


Superconducting 
Outsert Magnet 


Hybrid Magnet 
Cross Section 


Bore 


Resistive 
Insert Magnet 


Electrical ■ 
Supply 

Resistive 

Magnet 


Cooling 
Water Pipes 


Figure 2.6: Cutaway diagram of NHMFL 45T Hybrid magnet, from [21]. 


The superconducting “outsert” magnet wraps around the room temperature magnet. It 
is comprised of three connected subcoils in series with a design current of 10 kA. Each 
subcoil is made with of a niobium-tin alloy with a combined stored energy of nearly 100 
MJ. Many of the supported research areas for the hybrid magnet include: condensed matter 
physics, quantum hall effects, high temperature superconductivity, geochemistry, optical 
microscopy, resonant ultrasound spectroscopy and NMR studies, for example. The system 
was design for maximum user-friendliness, easy access to bore experiment area, and a 


19 


















10-year life with 2000 charge/discharge cycles for the “insert” or inner restive magnet [22]. 


For use within the laser accelerator system, static and continuous magnetic fields may not 
be needed because the time-scales of interaction are very short, typically < 10 ns. Pulsed- 
field systems are able to provide long duration pulses compared to the time-scales of inter¬ 
action for the laser accelerator are in effect “static." NHMFL Pulsed-Field Facility (PFF) 
at Los Alamos National Laboratory (LANL) offers researchers access to various capacitor 
driven magnets up to fields of 100 T (see Table 2.3) [19]. The fields can be maintained over 
a length of about one meter. 

The PFF system (see cutaway Fig. 2.7) combines a 40T outer magnet with a 60T field 
created by a capacitor-bank insert. The power system consists of a 1.43 GW generator 
combined with seven 64 MW power supply modules [23] driving each of the magnet’s 
seven coils. 

NHMFL Pulsed Field Facility is currently designing the next generation of non-destructive 
pulsed field magnets. The 100T magnet is a joint project between the National Science 
Foundation and the Department of Energy. Pulses up to 100T for periods of milliseconds 
will be available for users. Table 2.4 recaps the systems we have used in our examples. 

Table 2.3: Available Pulsed Magnets at PFF LANL. SP = Short Pulse. MP = Mid-Pulse, after [19] 


Cell No. 

Magnet 

Pulse Duration (ms) 

Bore (mm) 

1 

50 T SP 

25 

24 

2 

40 T MP 

400 

24 

2 

50 T MP 

400 

15 

3 

60 T SP 

25 

15 

4 

65 T SP 

25 

15 

5 

60 T SP 

25 

15 


Table 2.4: Magnet Specification Summary, after [19]. 


System 

Power 

Cooling Power 

Volume 

45T Hybrid 
Pulsed 

25 MW (continuous) 
1.43 GW @ 85 ms 

9 MW 

5 MW 

175 m 3 
35 m 3 
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Figure 2.7: Cutaway diagram of NFLMFL 100T pulsed magnet system, after [19] 
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CHAPTER 3: 
Equations of Motion 


3.1 Relativistic Lorentz Force 


We begin the laser accelerator analysis by considering whether a classical or quantum me¬ 
chanical treatment is required. The injector gives the electrons an energy of about 120 MeV 
with a one per-cent energy spread as input to the laser accelerator. This corresponds to a 
spread in electron wavenumbers of 

Ak e = « 6 ■ 10 12 m' 1 . (3.1) 

he 

So the corresponding uncertainty in position is 

Ax e > —« 8 • 10~ 14 m. (3.2) 

2A k 

The uncertainty in position is much smaller than the laser wavelength of 1 gm. The laser 
accelerator is clearly in a classical regime per Ehrenfest’s theorem [24]. 

Continuing the analysis with the relativistic Lorentz force law in the general case, an elec¬ 
tron moving through a magnetic field will spiral around the field lines. Without an applied 
electric field, the electron will not gain energy. The strong electric fields found in today’s 
high powered lasers can provide huge accelerating gradients. If a laser’s electric field can 
be coupled to an electron’s spiral motion, it can lead to significant energy gain. In compact 
covariant notation, the relativistic Lorentz force law [25] in SI units is 


dp ** dU» 
—— =m—— 
d T d x 


-eF^ v U v 


(3.3) 


where = (ymc. p) is the four-momentum, 7=1/y/l — (v/c) 2 is the Lorentz factor, m is 
the mass of the electron, U^ = y(c,v) is the covariant form of the electron’s four-velocity, 
r is proper time, e is the magnitude of the electron’s charge, and F^ v is the field-strength 
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tensor, 


pW 


f 0 -E x /c 

E x /c 0 

Ey/c B z 

\ E z /c —By 


Ey/c E z /c ^ 

—B z By 

0 -B x 

B x 0 j 


(3.4) 


We can transform between covariant and contravariant forms of 4-vectors in the usual way 
using the metric tensor 




(1 

0 

0 

0 ^ 

0 

-1 

0 

0 

0 

0 

-1 

0 

1° 

0 

0 

- 1 J 


(3.5) 


The left-hand side of Eq. (3.3) is a proper time derivative. In order to transform this into 
a lab-frame time derivative, we make use of the invariance of the spacetime interval. For 
some inertial frame K, the spacetime differential interval is given by 


ds 2 = c 2 dt 2 — dx 2 . 


(3.6) 


For another inertial frame K' . the interval is 


ds' 2 = c 2 dt' 2 — dx!^. 


(3.7) 


For sufficiently small spatial intervals, 


dx = \(t)dt. 


(3.8) 


Substituting Eq. (3.8) into (3.6) yields 

ds 2 = c 2 dt 2 { 1-/3 2 ) (3.9) 

where j3 = (J3 X , /3 V . /T) and /3 2 = v 2 /c 2 . Since v < c in all situations, ds 2 is always greater 
than zero in Eq. (3.9). For the given metric in Eq. (3.5), positive spacetime intervals are 
called timelike. Timelike intervals are ones where there exists a proper time frame K' such 
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that the spatial interval, dx', is zero. So we have from Eq. (3.7), 


ds' 2 = c 2 dt' 2 = c 2 dz 2 . 


(3.10) 


The spacetime interval is an invariant. Therefore, Eqs. (3.6) and (3.10) are equivalent. 
Combining the right hand sides of Eqs. (3.9) and (3.10) we now have, 

c 2 dT 2 = c 2 dt 2 {l-p 2 ). (3.11) 


We can now perform a change of variables in Eq. (3.3) using dt/dz = 7 from Eq. (3.11), 


d dt d d 
dz dz dt ^ dt 


(3.12) 


Therefore, Eq. (3.3) becomes 

^f =r ^ = . eF ^ Uv . (3.13) 

dz dt 

We will solve Eq. (3.13) subject to the laser field and an external uniform magnetic field 
where the appropriate forms for the electric and magnetic fields will be added. 


3.2 Equations of Motion 

Consider a circularly-polarized plane wave given by the following expressions for the laser 
electric and magnetic fields, E r and B r , with an external uniform magnetic field, B u , aligned 
along laser propagation axis,z, 


with 


E 0 (cos £, sin £, 0 ), 

(3.14) 

— (sin£,cos£, 0 ), 
c 

(3.15) 

( 0,0 ,B U ), 

(3.16) 

^ —— kfZ COyt. 

(3.17) 
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The plane-wave phase is given by £, the angular frequency of the laser light is given by co r , 
and the laser light wavenumber is given by k r . Substituting Eqs. (3.15) into (3.4) yields 




( 0 

-^-cosC 

yr sin £ 

° \ 

y 2 COS C, 

0 

-B u 

7? cos C, 

-'I" sinC 

B u 

0 

-^sinC 

V o 

COS C, 

f^sin C 

o / 


(3.18) 


for the field-stress tensor. Substituting Eq. (3.18) into (3.3) and running through compo¬ 
nents /i = 0 , 1 , 2 ,3 results in the following equations: 

[& cos £ - & sin £] (3.19) 

dt me 

4(rft) = - —(l-A)cosC-—ft (3.20) 

dt me m 

/] pf 1 

^(rJ8y)= — (l-/3,)sinC + — fix (3.21) 

at ■ me m 

d pF 

-T t (rPz) = -- [fix cos C — Ay sin £] (3.22) 

dt me 

where /3 X = v x /c, /3 V = v-y/c, and /3, = v z /c. Let us describe the geometry carefully (see Fig. 
3.1). Given initial electron energy y 0 we have the following definitions: 


Po = 
Pzp = 
Pxo = 
Pyo = 



cos a = 

P ± 0 cos (j ) 0 
P±o sin (j ) 0 


1 

7o 


77\/YZ- lcosa ’ 


P 0 sin a cos </> 0 
/3 o sin a sin (j ) 0 - 


1 

72 

1 

— i 

7> 



7r 


— 1 sin a cos (j) 0 , and 


1 sin a sin (j ) 0 


(3.23) 

(3.24) 

(3.25) 

(3.26) 


where a is injection angle of the electron measured from the z axis (optical and magnetic 
field axis), /3„ is the initial electron speed relative to the speed of light, j3 xo is the x compo- 
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nent of the initial electron speed relative to the speed of light (similarly for f5 yo and /3- 0 ), 
P± 0 is the perpendicular component of the initial electron velocity relative to the speed of 
light, and <j) 0 is the initial angle between and fi xo . 



y 


laser light 


Figure 3.1: Components of the laser accelerator interaction. The electric field vector is given by 
E r . Velocity vector (in units of c) is described by J3 with components Px’ Py’ and 
j3_. Transverse velocity is given by /3 ± . The injection angle is given by a. The angle 
between the rotating electric field and P± is defined as 6. The initial azimuthal angle, 
is defined as the angle between the x axis and J3 ± at t = 0. The external external 
magnetic field is given by B u . 


Notice that the right-hand side of Eqs. (3.19) and (3.22) are equivalent. Equating the 
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left-hand sides to each other and integrating gives a constant of motion, 


dy 

dt 



7(1 ~Pz) = e 


(3.27) 


where e is a constant for all time [1], Both y and /E are dynamical variables. As /increases, 
the factor 1 — j3 z must decrease to keep £ constant. We can express /T 2 = /3 2 + /3 2 in terms 
of £ as follows. From the definition of the Lorentz factor we have, 


7 = 




l-P?-ft 


/5i = ^7-i-fft 2 ]- 


(3.28) 


Using Eq. (3.27) for /3. 2 results in 


/5± — 7 [2ye — £ 2 — l] ■ (3.29) 

From Eqs. (3.27) and (3.29), as the electron gains energy from the laser accelerator, y 
increases. Therefore, /3- will increase while f3_ decreases. The coupling to electric field is 
governed by /3j_. 


Now we will attempt to uncover other constants of the motion by extracting a full time 
derivative. Consider the following relations, 


d 

dt 


sin£ 


— sin(k,z — (O r t) 

cos (k r z — cod) — ( k r z — cod) 
dt 

cos (k r z — cod) (. k r z — (Of) 
—(Of( 1 — /E) cos (k r z — (Oft) 
~(0f( l-/3 z )cosC 


(3.30) 


and similarly, 


d 

dt 


cos £ = Q) r ( 1 — fi z ) sin C,. 


(3.31) 
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Furthermore, applying the chain rule on d/dt{y[3 X sin Q and substituting Eq. (3.30) yields 


j f (YP.x sin C) = sin C j^YPx) + Ypxj f (sin Q 

= sinC^W-rWl-^cosC, 


(3.32) 


and also using Eq. (3.31), 



^(7Av COS O = COS C J t (Ypy) + rPyCOri 1 ~ Pz) Sm £. 

(3.33) 

A second constant of motion can be determined by solving Eqs. (3.20) and (3.21) for f3 x 

and /3 V : 

«■= 

(3.34) 



(3.35) 


Substituting these equations into Eq. (3.19) and making use of the relations in Eqs. (3.30- 
3.33) results in 

T, = -^,Jl {Y ^ COsQ ~ + J^rd - ftHtftsinC-ftcosO. 

(3.36) 

Using Eq. (3.19) to solve for (/3 v sin£ — p x cosQ in terms of dy/dt,Eq. (3.36) becomes 


dy E () d . _ E () d . _ . m _ . dy 

-T- + — — (r/3yCosC) + — — (y&sm£)- — y(l-/3 z )w r — = 0. 


(it cB„ (it 


c5„ dt 


eB,, 


dt 


(3.37) 


Extracting an overall time derivative gives 


d 

dt 


E (> E () . m 

7+ -^-ypy cos C + —rftcSinC - — ea»,7 
cB u 


- 0 . 


(3.38) 


Factoring y and making use of the constant e, defined in Eq. (3.27), yields 

= 0 


d 

dt 


F F /17 

7(1 + -^/5yCosC + - 777 - /3 a sin £ — — y 0 (l -p zo )(O r 


(3.39) 
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where y () is electron’s initial Lorentz factor value at / = 0, and e’s value can be set from 
initial conditions such that £ = y 0 { 1 — /3- 0 ). Multiplying both sides by eB u /my„, we can 
integrate exactly, 


7 


pH pH 

— ~M 1 - Pzo) + Z~(Px sin £ + Py C0S 0 

y 0 m y 0 mc 


= constant 


In Eq. (3.40) we can define the following constants, 

eB u 


(O r 


y 0 m 

eE (> 


(Oe = 

y 0 mc 

Am = co c — o) r ( 1 — p zo ) 


(3.40) 


(3.41) 

(3.42) 

(3.43) 


where (O c is the relativistic cyclotron frequency, Aft) represents the difference between the 
relativistic cyclotron frequency and the Doppler shifted light frequency [1], and (Oe is elec¬ 
tric field analogue to the relativistic cyclotron frequency. With these substitutions, Eq. 
(3.40) becomes 

y[Aft) + (Oe (fix sin £+P y cos Q] — constant = 77 . (3.44) 


Two other constants of the motion can be developed from the transverse equations of mo¬ 
tion, Eqs. (3.20) and (3.21). The same procedure of eliminating (1 — /3 z )cos£ in favor of 
a time derivative can be applied on Eqs. (3.20) and (3.21). For example, substituting Eq. 
(3.30) into (3.20) yields, 


d 

dt 


( 7 / 3,0 


eE 0 

me 



eBu 

m 


Pr 


(3.45) 


We can factor out a time derivative and integrate the expression exactly giving us 


rPx- 


eE a . eB u 

-sin Q -\ - y = constant = a. 

me (O r me 


(3.46) 
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Similarly, using Eq. (3.31) in (3.21), we arrive at 


y/3 v -— cos C, - -x — constant = b. (3.47) 

mc(O r me 

Table 3.1 summarizes the equations of motion and constants of motion where the electron 
phase is £ = k r z— co r t. Our goal is to solve these equations for the electron’s position and 


Table 3.1: Equations and Constants of the Motion 


Eq. No. 

Differential Equation 

Eq. No. 

Constant of the Motion 

3.19 

Tt = — me [Ac COS C — Ay sin C] 

3.27 

£ = 7(1 ~Pz) 

3.20 

UrM = -^cosC(l-/3 z )-^j8 y 

3.44 

rj = Y[A(0 + (0 E (PxSin£ + PyCosQ] 

3.21 

&m= ^sinC(l-/3 z ) + ^/3, 

3.46 

a = rP* 

3.22 

|(rft) = ^ [ft cos f-ft, sin f] 

3.47 

b = rPy mX COS £ mc X 


energy and to determine the optimal conditions for significant energy gain. 


3.3 Integral Representation 

In order to solve Eqs. (3.19) through (3.22) and arrive at an expression that relates distance 
traveled in the interaction region to energy gained, it will be beneficial to recast some of 
the equations in Table 3.1. Consider, 


_l ■ E r = (/3 x ,/3 y ,0) -^(cosC, —sinC,0) 
= E 0 \fi x cos £ — /3 y sin Q 
= cos 0, 
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where P± = \J PI + /3 2 and 6(t) is the angle between the electric field vector, E r and P ± 
such that 0 < 9{t) < n. Therefore, 


[/3 x cos£ — p y sin Q = /3^cos0. (3.48) 

Also consider the following cross-product, 

p± X E r = (p x ,py,0) x Eo( cos£,-sin£,0) 

= -E 0 [p x sinC + PyCosQz. (3.49) 

Squaring both sides we have and using the magnitude definition of the cross-product we 
have, 

(0±x E r ) • Q5j_ X E r ) = E 2 [P x cos £ + P y sin C] 2 = P±Eg sin 2 0. (3.50) 

Therefore, 

[Px sin C + Py cos C] = P± sin 0. (3.51) 


With Eqs. (3.48) and (3.51) in mind, let’s return to Eqs. (3.19) and (3.44). We now have 


dy 

dt 


o n 

-Pj^COS 0 

me 

y{Aco + (OEp^sinO). 


(3.52) 

(3.53) 


Note that dt — ydz/c(y — £) from Eq. (3.27). Substituting Eqs. (3.29) and (3.42) into Eq. 
(3.52), we now have, 


dy = 


Yo(Q e [2yg-£ 2 - l] 1 / 2 
c (r-e) 


cos 9 dz. 


(3.54) 


We can eliminate the cos 9 factor by solving Eq. (3.53) for sin 9 and noting that sin 2 9 + 
cos 2 0 = 1, we then have, 


cos 0 


± 


ft)|(2ye-e 2 - 

®l(2ye 


1) - (yAco-p) 
-e 2 - 1) 


21 1/2 


(3.55) 
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Substituting the above equation in Eq. (3.54) we arrive at 


dy= ±Yo 


[coj(2y£ - e 2 - 1) - (yAco - rj) 2 ] 1/2 dz. 


c(y~£) 


(3.56) 


In utilizing the constants of the motion e and 77 , we have arrived at an equation that does 
not explicitly depend on the angle 0 between j3j^ and E r . Integrating both sides of Eq. 
(3.56) gives us 


y - £ 


r 0 [ft)l(2y / e-e 2 - 1) - (y'Aco-ri) 2 


1/2 


dy 1 



(3.57) 


where the negative solution was dropped since we are interested in energy gain. The left- 
hand side is a family of curves known as curtate cycloids [1], while the right-hand side 
yields the distance traveled along interaction region. 


3.4 Curtate Cycloids 

We can motivate a geometric understanding of Eq. (3.57) by considering the mathematics 
of cycloids. A cycloid is the curve formed by tracing a point either inside, outside, or 
on the rim of a circle of radius q rolling along a straight line [26]. The general equation 
characterizing the curve is written in parametric form: 

X = tqyr- k sint/r 

Y — — fccost/r (3.58) 

where q is the radius of the circle, K is the distance from the center of the circle from where 
the curve is traced, yr is the parametric angle. There are three types of cyloids, as shown in 
Fig. 3.2. Standard cycloids have K — curtate cycloids have K < tq, and prolate cycloids 
have K > <q. 

Consider a curtate cycloid with K < <q. If we take differentials of Eq. (3.58) then 

dX =(^ - Kcosy/)dy/ 

dY =k sint/r dy/ (3.59) 


33 



Standard Cycloid 



Curtate Cycloid 



Standard Cycloid 



Figure 3.2: Diagram describing different types of cycloids: standard cycloid, curtate, and prolate. 

The radius of the wheel is E,. The distance along the radius where curve is drawn is 
given by k, from [26]. 


Now consider the product Y dY, 

Y dY = — fccos \j/)(K sin \j/d\j/) — Ksiny/dX. (3.60) 


Therefore, 

_ YdY 

dX = — : -. (3.61) 

fcsin y/ 

We can eliminate the sin y/ term in Eq. (3.61) with the standard Pythagorean trigonometric 
identity and Eq. 3.58 resulting in 


dX 


YdY 

\J-Y 2 + 2E,Y + (k 2 - E, 2 ) 


(3.62) 
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Notice this expression matches the form of Eq. (3.57). Integrating both sides of Eq. (3.62) 
in order to get an explicit expression for x as a function of y, we have 




— y' 2 + 2<^y / + k 2 — <§ 2 — E, sin 1 


izi: 

K 


(3.63) 


We have to be careful when employing the limits of integration since the square-root term 
must be real and arcsine function is valid only when the argument is less than or equal to 
111. The arcsine requirement provides the following constraint on Eq. (3.63): 


K 


< 1 . 


(3.64) 


Squaring both sides of Eq. (3.64), we have 


-Y 2 +2%Y + k 2 -% 2 >0. 


(3.65) 


The radicand of the square-root term in Eq. (3.63) matches Eq. (3.65). Therefore, only one 
constraint on the value of Y exists. Solving the inequality gives the following constraint on 

F: 

-k<Y <% + k. (3.66) 

Let us consider a simple example where — 3 and K — 1 in Eq. (3.58). Therefore, the 
minimum and maximum values for Y are 2 and 4, respectively per Eq. (3.66). The explicit 
function for X using Eq. (3.63) is 

X(Y) = - \/-F 2 + 6F-8-3suW 
for 2 < Y <4. In order to find Y(X) we would need to 


1 (3-F) + ^ (3.67) 

invert Eq. (3.67) numerically and 
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use periodic boundary conditions such that minima and maxima of Y occur at 

3fmax = (2 tn-\- 1)7T<§ 

Y max = ^ K 

2fmin = 2mnE> 

Fmi„ = %-K (3.68) 


form = 0,1,2,3... [26]. 

In order to make the appropriate comparison to Eq (3.62) we must manipulate Eq. (3.57) a 
bit. Let y = y— £ and factor Act) from the radicand. Therefore, Eq. (3.57) becomes 



y'dy 



+ 2 


a>££+2riAa)-2A(0 2 £ 

Ad) 2 


y+ 


co|e 2 -(o| 


r] 2 +2riAco£~£ 2 Aco 2 
A co 2 


A CD—Z. 


(3.69) 


We can now make the following identifications. The y axis is proportional to the energy of 
electron and the abscissa is proportional to the z distance traveled by the electron. Relating 
this to the curtate cycloid, the “radius” of the wheel is related to the amplitude of energy 
oscillations as the electron moves down the laser accelerator. The radius is given by 


cu|£ + 2 t 7 Aco - 2A are 
Acu 2 


(3.70) 


and the point from the center of the wheel where the curtate cycloid curve is drawn is given 
by 


K = 


cole 2 -co£-ri 2 + 2r\ Acoe - e 2 Aco 2 


A ft) 2 


+ 


ft)|£ 2 + 4?7“Aft) 2 + 4Acu 4 £- + 4co|£r7 Acu 


4ft)-Acu 2 £ 2 - 8rjAco^£ 


A ft) 4 


(3.71) 


An explicit equation for distance as a function of energy gained can be immediately identi¬ 
fied through Eq. (3.63) with appropriate substitutions for and K given by Eqs. (3.70) and 
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(3.71) with the constraints on y such that <§ — k <y— £ < ^ + K with the form 


z 


— y /2 + 2E, y' + k 2 — t, 2 — sin 


r 


K 



(3.72) 


3.5 Resonance Condition 


The resonance condition is defined as 


Aft) = 0 

Aft>= —-ft),.(l-/3-o) =0 ( 3 . 73 ) 

77» 


This is the case where the electron’s relativistic cyclotron frequency is equal to the Doppler 
shifted light frequency. When these two frequencies are balanced, the electron experiences 
a sustained acceleration. Depending on initial conditions, the electron may initially gain or 
lose energy. Furthermore, if resonance condition satisfied, then it is maintained throughout 
the interaction. Since e is a constant of the motion, we can rewrite Eq. (3.73) as 


eBu E 

-ft), - = 0. 

ym y 


(3.74) 


Multiplying by yon both sides of Eq. (3.74), we have 


eBu 

-ft), £ = 0. (3.75) 

m 


Everything on the left-hand-side of Eq. (3.75) is a constant, therefore Aft) is zero for all 
time. 


There are infinitely many combinations of y 0 , B u , A,-, and /3 - 0 in Eq. (3.43) that would drive 
A ft) to be zero. We can use the resonance condition as a constraint to help us sample the 
parameter space. 
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3.6 Electron Phase 


Energy gained or lost by electrons in the laser accelerator is given by Eq. (3.52). Recall 
that 0, a dynamical variable, is the angle between J3j_ and E r as shown in Fig. 3.1. If 
0 — K, then the right-hand side of Eq. (3.52) is positive since cos(^) = — 1. Therefore, y 
will grow with time. If 0 — n/2, the right-hand side is zero, y will be at an extremum. The 
seemingly worst case for a laser accelerator system is when 0 = 0, then the right-hand side 
of Eq. (3.52) will have the greatest negative value resulting in electrons losing energy. The 
initial angle between J3 ± and E r at the start of the interaction is defined to be 0 o . 

The phase for each electron given by £ = k,-z — (O r t will have a value an initial value of £ 0 = 
k r z 0 at the start of the interaction since t = 0. Experimentally, the laser accelerator would 
be accepting relativistic electrons from an injector system. There will be a distribution of 
electrons spread along z as shown in Fig. 3.3. This spread in is equivalent to 0 < Q 0 < 2 n 
over a wavelength of light. So we expect some electrons to gain energy; some electrons 
to lose energy. But as we will see in Chapter 5, having an unfavorable Q„ may not be so 
detrimental for achieving good electron energy gain over the interaction length. 


Electric 



Figure 3.3: Diagram of initial electron distribution in the longitudinal direction over two wave¬ 
lengths of light. In reality there would be millions of electrons in the distribution, af¬ 
ter [27]. 
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3.7 Maximum Energy Gain 


The ideal case for maximum electron acceleration and energy gain is when the resonance 
condition is met or when Aft) = 0 and with 0 o — K. With this substitution, Eq. (3.57) 
becomes, 


T 


dy' = 


(3.76) 


7o 


ft)|( 2 y / e-e 2 - 1 ) -r\ : 


c 


Integrating both sides and solving for z we have, 


z = 


Yo 3co|e 2 [ 

For large values of y, z y 3//2 [1]. 


2 co|ey / - G)|e 2 - co| - T 7 2 ( (Dfey' + ft)|e 2 + co| + 77 2 - 3ft)|e 2 


(3.77) 


As mentioned before, there are infinitely many combinations of y 0 , B u , X r , and fi z0 for 
resonance to be met. Practically, an experiment can be setup by choosing the laser and an 
external magnetic field leaving y 0 and fi z0 as free variables in order to find resonance. At 
resonance using Eq. (3.41) we have, 


A ft) 


eB u 

my 0 


®r(l fizo ) 


0 


so, 

pB 

e -l = co r (l-p zo ) 

my (> 


Using Eq. (3.27) we have 


eB u 

6) r m 


(3.78) 


With Eqs. (3.24) and (3.78), noting that f5 z0 = (y 0 — e)/y 0 again from Eq. (3.27), we can 
solve for the injection angle 


a = cos 


Yo-e 

Vvh 1 


= cos 


To 


eBj, 

CO,m 



(3.79) 


In order for a to be a real-valued angle, there’s a requirement on the inverse-cosine function 
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of Eq. (3.79) such that the argument must be less than or equal to one, 


Yo-e 

Vyo- 1 


< l. 


(3.80) 


Squaring both sides and completing the inequality for y () gives us a minimum to achieve 


resonance, 


7min 


e 2 +1 
2e 



(3.81) 
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CHAPTER 4: 
Numerical Methodology 


4.1 Finite-difference Method 


Both Eqs. (3.72) and (3.77) represent analytical solutions for the off-resonance and reso¬ 
nance cases, of distance along z as a function of /respectively; but they are quite unwieldy. 
Finite difference numerical integration methods provide a much cleaner solution. 

We can rewrite Eqs. (3.19) through (3.22) in terms of four-momentum components. Re¬ 
calling that pP = my(c,v x ,v y ,v z ) we have, 


dp 0 

dt 

dp 1 
dt 

dp 1 
dt 

dp 3 
dt 


eE ° r„i 


[p cos£ — p~ sin £] 


-eE 0 1 


r eB u 2 

cos C- p 

ym 


eE 0 1--E- sinC + — p 
\ ymc I ym 

— eE d r/cosC-^sinCl 

ymc L J 


(4.1) 

(4.2) 

(4.3) 

(4.4) 


where £ = k r z — co r t, p° = /me, p 1 = ymfi x c, p 2 = ymj3 y c, and /; 3 = ymfi z c. Turning the 
first equation into a finite-difference equation, we have 


E old + 


- [p x cos £ - sin Q At 

7 mc ■ J old 


The “old” and “new” labels corresponds to previous time step values and updated time step 
values, respectively, of the variable in question. We can then determine the value of y new 


through 


/new 


P new 
me 


(4.6) 
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Then using Eq. (3.27) 


P 


new 

z 


"/new £ 
/new 


since £ — y(l — /3-) is a constant of the motion. The position in z is then 


(4.7) 


Znew ~ Zold + ft?™ C At. 


(4.8) 


Similarly for the x component, 


P 


l 

new 


7’old + 


-eE 0 



cos £ 


eBu ? 1 
- p At 

Y m J old 


Then 


Therefore the x-position is 


A „ew 


P new 
TtiewinC 


-Tiew ~ -^old T P x ^^ • 


(4.9) 


(4.10) 

(4.11) 


and similarly for the y component. Since we are updating the new positions on the updated 
velocities, this is known as Euler-Cromer method [28]. 


Numerical integration using a 4th order Runga-Kutta (RK4) method was also performed. 
The percent difference in electron energy gain between RK4 and Euler-Cromer methods is 
only 0.03 percent. Since the percent difference is negligible, Euler-Cromer methods were 
employed in the dissertation due to their speed and simplicity. 


4.2 Initial Conditions on Resonance (Aco = 0) 

The task in this section is to numerically solve the equations in Table 3.1 subject to the 
electrons entering the interaction region with the resonance condition satisfied such that 
A co = 0. We ultimately want to know the energy gained by the electron and its trajectory. 

Assume that the initial electron energy and external magnetic field are independently cho¬ 
sen parameters. From Eq. (3.79), we have an equation for the injection angle, a — 
a(y 0 ,B u ), as a function of two free variables. The following result shown in Fig. 4.1 is 
a surface of required injection angles for resonance. 
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Figure 4.1: Shown is the the required injection angle for resonance as a function of initial electron 
energy, y c , and external magnetic field, B u , for a given ft),- 


If we select a magnetic field strength of 45 T, then we can slice through a constant external 
magnetic field in Fig. 4.1 resulting in Fig. 4.2. We now have an optimal injection angle 
for various initial electron energy values where the minimum value of y 0 in Fig. (4.2) is 
consistent with Eq. (3.81). 
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Required injection angle for B u = 45 T 



Figure 4.2: Shown is the allowed ordered pair combination of initial electron energy, y, and initial 
injection angle, a for a given magnetic field to guarantee resonance for a given co r . 


Consider the electron at the start of the laser accelerator system at time, t — 0. Using Eqs. 
(3.15) and (3.17), the electric field, E r , is then oriented along the jc axis or E r = E a x since 
£(f = 0) = k r z — (O r t — 0. From Eq. (3.52) it would seem that the greatest energy gain 
would require an anti-alignment between E r and /3 ± such that 0 — K initially. The other 
spherical coordinate angle, </> 0 , which is the angle between /3 ± and f$ xo , is then also equal 
to n. In general, <f) 0 may take on any value between 0 and n. So, with a and 0 O established, 
the initial conditions for position and velocity can be calculated using Eqs. (3.23) through 
(3.26). 

4.3 Electron Trajectory 

For plane-wave propagation (no diffraction) with 0 — % for maximum acceleration, both 
Eqs. (3.52) and (3.56) are equivalent with A ft) = 77 = 0. It is straightforward to numerically 
integrate and update y. Then for each timestep using the constant of the motion e, /3- can 
be updated and the position z calculated. The values for /3 t and /3 y can also be evaluated by 
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making use of constants a and b where their values are given by initial conditions. From, 
Eqs. (3.46) and (3.47), /3 V and f5 y are 


Px 

Py 


eE 0 . cB u 

a H-sin Q - y 

me (O r me 

, eE 0 y. eB u 

b H-cos Q H-x 

me (O r me 


(4.12) 

(4.13) 


We can determine the number of times an electron will spiral around the external magnet’s 
field lines by integrating Eq. (3.41) (and making the substitution of y{t) for y Q ) over the 
time it takes the electron to clear the interaction region. Stronger magnetic fields results 
in tighter and more rotations. Trajectories from numerical integration will be shown in 
Chapter 5. 


4.4 Synchrotron Radiation 


An estimate of power loss through synchrotron radiation can be made if we consider the 
simplified case with the laser source switched off. With only the external magnetic field 
and no laser field present, E 0 = 0, Eqs. (3.19) through (3.22) reduces to 


7(0 — 7(0) = 7)7 constant 

dv x _ _eBu_ 
dt my y 
dv y eB u 

= - v X 

at my 


v z(0 — V z (0) = v zo , constant 


(4.14) 

(4.15) 

(4.16) 

(4.17) 


Eqs. (4.15) and (4.16) can be solved by letting p = v x + iv y and (O c — eB u /my. Therefore, 
p = i(O c p. Solving this with initial conditions yields p(t) = vj_ 0 e4 COff + 5 ) where v± 0 = 
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v xo + v yo an d <5 = tan 1 ( v ya / v xo ). Integrating once more yields the following 


x(t) = 

y 0) = 

z(t) = 


—^sin(co c t + d) +x 0 — —— sin(<5) 
CO c CO c 

-—^-cos(co c t + 8) + y 0 + —— cos(<5) 

(O c CO c 

Zo "I” V zo t. 


(4.18) 


The radius of the orbit is then given by the cyclotron radius for an electron in a static 
magnetic field. 


r= \/(,x(t)-x 0 ) 2 + (y(t)-y 0 ) 2 
^ ymv±o 
eB,, 


(4.19) 


We can calculate the synchrotron power loss using the relativistic Larmor formula given 
by [29], 


P = 


A i 0 e 2 fc 
6 n 


p 1 - Pxfi 


(4.20) 


Taking one and two time derivatives of Eqs. (4.18) and substituting into Eq. (4.20) yields, 




6 n 




P = 




6 71 


We arrive at [29], 


floe Yc 202 
p - 6 k ( ° cFa -°- 


(4.21) 


Plugging in nominal values for an external magnetic field of 45 T and 300 T, we can expect 
synchrotron energy losses to be on the order of several electron-volts (eV) to 5 keV, respec¬ 
tively. These figures represent an upper-bound estimate energy losses. Once the laser field 
is present, as will be shown later, y increases however J3j_ is small and its value decreases. 
For energy gain of approximately one GeV, power loss through synchrotron radiation is 
negligible. 
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4.5 Numerical Integration 

The numerical integration procedure for off-resonance or resonant analysis begins by se¬ 
lecting global constants such as laser wavelength, laser power, estimated spot size, external 
magnetic field strength, and time steps. If resonance analysis is required, the angle for res¬ 
onance based on the global constants is calculated. We compute initial momentum values 
in preparation for numerical integration. Based on the initial values, we perform Euler- 
Cromer integration to determine “new” values for y and the momenta. We use these new 
values to calculate position and the overall trajectory of the electron. We also calculate 
power loss through synchrotron radiation and perform diagnostic checks of the code for 
each time step. Finally, we keep track of the values by writing the data to a file for post¬ 
processing in MatLab. 
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CHAPTER 5: 

Numerical Analysis Results 


In this chapter, we will use the numerical integration method previously described to obtain 
simulation results on the laser accelerator performance for a variety of initial conditions 
including: on resonance, different values of the initial phase 0 o , the angle between j3 ± and 
E r , off resonance, and real field effects. 


5.1 On-resonance with 0 o = n 

Recall that resonance is defined as A 0 ) = 0 and occurs when the electron “feels” the Doppler 
shifted circularly-polarized laser light oscillating at the same frequency as the electron’s 
spiraling motion in the uniform magnetic field. In the simulation, we assume that the 
external magnetic field B„ and laser source parameters (laser power and wavelength) are 
fixed. The resonance condition is then a function of the injection angle a and injection 
energy y () as shown in Fig. 4.2. 

Let us assume the following for the simulation of a single electron . The laser source will 
be like the Texas Petawatt Laser (TPL) system. The laser Rayleigh length, Z,, is assumed 
to be longer than interaction region of one meter in order to eliminate diffraction effects. 
Gaussian beam analysis will be considered later in the chapter. The optical mode waist is 
about 1 mm throughout the interaction region. The external magnetic field is assumed to 
be static, uniform field throughout the interaction region. Strong magnetic fields support 
electron to light coupling for energy gain. Let us also assume 9 0 = K so that /3j_ is anti¬ 
parallel to the electric field of the light as shown in Fig. 3.1. 

For the single electron simulation, Table 5.1 summarizes the input parameters. 

Performing the numerical integration as described in the previous chapter, we obtain the 
results shown in Fig. 5.1. The results are compared to the analytic solution given in Eq. 
(3.77). The percent difference is about 0.2 percent for the final energy of the electron after 
a one meter long interaction. We compared multiple simulation results to the analytic pre¬ 
diction in regimes where the analytic formation remained valid. In all cases, the simulation 
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Table 5.1: Input Parameters for Resonance Laser Accelerator Simulation 


Variable 

Symbol 

Value 

Initial Electron Lorentz Factor 

Yo 

119 

Initial Electron Energy 

<§o 

61 MeV 

Laser Peak Power 

p 

1 PW 

Laser Wavelength 

A,- 

1 gm 

External Magnetic Field 

B u 

45 T 

First Constant of Motion 

e 

4.2 x 10 -3 

Initial Azimuthal Angle 

0 o 

K radians 

Injection Angle 

a 

20.8 /i radians 

Time Step 

At 

50 fs 


always closely matched the analytic prediction. 



z (m) 


Figure 5.1: Comparison of analytic solution to numerical simulation results of an electron in the 
laser accelerator. Energy gain versus distance is plotted. The black curve describes the 
analytic solution and the red curve describes the numerical solution. The difference is 
about 0.2 percent for the final energy gain, which is approximately 1.17 GeV. 
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Now that we have established good agreement between analytic and numerical solutions 
we can now explore how other parameters such as electron radial position and velocity 
evolution over the interaction region. Fig. 5.2 shows simulation results for two electrons 
at two different external magnetic field values, 45 T (blue) and 100 T (red) with injection 
angles a of 20.8 /irad and 9.3 mrad respectively for resonance as per Eq. (3.79). The other 
parameters are the same as in Table ??. 



Figure 5.2: Electron simulation results for the parameters given in Table 5.1. The blue curve give 
results for B„ = 45 T while the red curve give results for B„ = 100 T. The top left 
graph plots electron energy versus distance traveled. The top right graph plots the radial 
distance of the electron from the optical axis. The bottom left graph plots f5± versus 
distance traveled. The bottom right graph plots 1 — versus distance traveled. 


The top left curve shows that the electron’s energy increases to nearly 1.2 GeV for 45 T 
and 1.6 GeV for 100 T in one meter. The electron’s orbit spirals out to a radius of about 
0.17 mm for 45 T and 0.15 mm for 100 T as shown in the top right graph. The value of 
1 — /3- is plotted in the bottom right curve since /T is nearly unity for the entire length of the 
interaction region. Note that f$ z increases while /3 decreases during most of the interaction. 
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In the bottom left curve, we can see that /3_ has a peak then slowly decreases for the 45 T 
case while the 100 T does not have a peak. Let us try to understand why /3j_ has a peak for 
the 45 T case while no peak exists for the 100 T case. It is not immediately clear from Eq. 
(3.29) alone that /3j_ should have a maximum. Calculating its derivative we have 


_ d_ f 1 
dz dz |_y(z) 

-I A 


2y(z)e-e 2 -l 


[2y(z)e-e 2 -l] 


2 ,1-1/2 i „2 ^n-\dy(z) 


~\r(z)~ ^ f(z) 

Setting Eq. (5.1) equal to zero we have an extremum when 

* e 2 +1 
7 = — —• 


m e-e 2 - ir -^ 


(5.1) 


(5.2) 


Notice that y* = 2y m i n from Eq. 3.81. Calculating y* for 45 T and 100 T cases we have 


y*(45 T) « 238 
y*(100T) « 107 


The value for y*(100 T) explains why there is no peak for the 100 T case in Fig. 5.2 since 
the simulation began with an initial y 0 = 119 T. Substituting Eq. (5.2) back into Eq. (3.29) 
yields 

fl_L,max = ; 9 =, (5.3) 

Ve 2 + 1 

which agrees with the simulation value for /3_ >max . In Fig. 5.3, we have again plotted the 
electron trajectories for two the cases: 45 T field in the blue curve and 100 T field in the 
red curve. Notice as y increases along the z axis it causes the radius of spiraling motion 
around the magnetic field line to increase. There are more numerous spirals for the stronger 
field. The number of spirals gives an approximate number of wavelengths of light that pass 
over the electron near resonance. High powered lasers such as the Texas Petawatt Laser 
produces Gaussian pulses. The Fourier transform of these pulses would also have a spread 
of frequencies in the frequency domain. As shown in Fig. 5.3, only about ten wavelengths 
of light pass over the electron. To first order we do not need to be concerned with the effects 
to the frequency changes associated with the pulse shape. 
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Figure 5.3: Electron trajectories in the laser accelerator for two different external magnetic Field 
CasesLaser accelerator for the parameters given in Table 5.1 in the blue curve. The 100 
T field trajectory is shown in the red curve for comparison. 


Table 5.2 summarizes the simulation results for both cases. 

Table 5.2: Results for Single Electron Simulation of Laser Accelerator at 45 T and 100 T 


Variable 

Symbol 

Value at 45 T 

Value at 100 T 

Final Lorentz Factor 

7/ 

2282 

3109 

Final Energy 


1.17 GeV 

1.59 GeV 

Final Electron radial position 

r 

0.16 mm 

0.13 mm 


5.2 Electron Phase Performance 

A realistic electron bunch is thousands of times longer than a wavelength of light. As 
mentioned in § 3.6, there will be millions of electrons in each optical wavelength that are 
randomly distributed in the microscopic longitudinal direction. Recall that £ = k r z — (0,t 
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and at time t — 0, £,,, = k r z 0 — 2 nz 0 /h- Their initial £ 0 ’s will therefore span a 2k range 
over a distance where z 0 changes by a wavelength of light, X r . We expect some electrons 
to gain energy and some electrons to lose energy depending on £ 0 . 

We select the coordinate system so that electron is initially at the origin and the laser electric 
field is initially along the x-axis. In order to sample different microscopic values of £ 0 , we 
consider three representative values of 9„ = 0. k/I.k in the simulations. Recall from Eq. 
(3.52), that these values are anticipated to cause the electron to initially lose or gain energy. 

Consider the following simulation, shown in Fig. 5.4 with y 0 — 130 (above the minimum 
energy required for resonance), B M = 45 T, and on resonance injection with a res « 2.4 mrad. 

In Fig. 5.4, we examine the evolution of energy versus distance traveled for the first 15mm. 
The middle graph shows /3j_ versus distance. The bottom graph shows the evolution of 
0. The colored curves shows the evolution of the respective parameter for three different 
initial 9 0 . The black curves shows 9 () = n, the green curves show 9 0 — n/ 2, red curves 
show 9„ = 0. 


In describing the energy plots we must recall Eq. (3.52) 


d 7 eE o a n 

— — -P^COS0. 

dt me 


Notice the factor of cos 0 on the right-hand side of the energy equation. For the 9 0 — n 
case shown by the black curve, the energy initially increases since the right-hand side of 
Eq. (3.52) is positive. For the the 9„ — k/2 case shown by the green curve, the energy is 
initially constant since the right-hand side of Eq. (3.52) is zero. For the 9 0 — 0 case shown 
by the red curve, the energy initially decreases since the right-hand side of Eq. (3.52) is 
negative. 


Fet us now look at the perpendicular velocity graphs for /3j_. For both the 9 0 = n and 
9 a = k/2 cases, J3j_ goes up initially, as shown by the black and green curves, consistent 
with the electron’s increasing energy from the top graph. For the 9 0 = 0 case shown by 
the red curve, J3j_ decreases in value until it reaches zero. When goes to zero, the angle 
0, shown in the bottom graph, changes from 0 to K quickly. Then f3 and the energy both 
begin to grow. This is a surprising result. It was not expected that an electron that starts out 
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— e 0 =jt — e 0 =7c/2 — e 0 =o 




o' -- 

0 5 10 15 

z (mm) 

Figure 5.4: Simulation results for on-resonance but different initial d 0 for y 0 = 130, B„ = 45. The 
top graph shows energy versus distance traveled for the first 15mm. The middle graph 
shows versus distance. The bottom graph shows the evolution of 6. The black curves 
shows 0 o = n. The green curves show d a = n/2. Red curves show 0 o = 0. 


moving parallel to the light’s electric field would eventually gain energy. 

Let us consider a similar analysis of 0„ but with lower initial electron energy of y a — 125 
and a res « 1.8 mrad. These results are shown in Fig. 5 . 5 . The electron’s perpendicular 
velocity / 3 _ starts out with a smaller value than the y Q — 130 case since the required injec¬ 
tion angle is smaller. So for the 9 a — 0 case, / 3 _ goes to zero sooner along the z axis. This 
causes 6 to change to K sooner than the y a — 130 case. We again see that an electron that 
starts the interaction with an unfavorable 6 0 , eventually changes to a favorable 0. Thus 
they accelerate to high energies along with electrons that initially started with favorable 
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Figure 5.5: Simulation results for on-resonance but different initial 0 o for y 0 = 125, B„ = 45. The 
top graph shows energy versus distance traveled for the first 15mm. The middle graph 
shows versus distance. The bottom graph shows the evolution of 6. The black curves 
shows 0 o = 71. The green curves show 6 0 = n/2. Red curves show 0 o = 0. 

9 0 values. Again, this is a surprising and encouraging result. Initially, it was thought that 
electrons that did not have a favorable 0 o would not gain energy. 

Let us again consider a similar analysis of initial 0 o but with an even lower initial electron 
energy of y 0 — 120 and a res ~ lmrad. These results are shown in Fig. 5.6. With y a — 120 
case, /3j_ starts out with an even smaller value than the previous two cases since the required 
injection angle is smaller. Therefore, j3± goes to zero sooner which causes 0 o to change 
quickly n sooner. 

Let us again consider a similar analysis of initial 0 o but an initial electron energy of y a — 119 
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Figure 5.6: Simulation results for on-resonance but different initial 9 0 for y 0 = 120, B„ = 45. The 
top graph shows energy versus distance traveled for the first 15mm. The middle graph 
shows j3|_ versus distance. The bottom graph shows the evolution of 6. The black curves 
shows 9 () = 71. The green curves show 6,, = n/2. Red curves show 0 o = 0. 


and a res « 20 /irad. These results are shown in Fig. 5.7. We are now quite familiar with 
the story. In the middle graphs, /3_ starts with a much smaller value than the previous three 
initial energy cases. For the 0 o — 0 case, /3± evolves to zero even sooner which quickly 
leads to 9 — 7T. This result encourages us to inject at the minimum energy possible where 
the resonance condition will remain satisfied, given in Eq. (3.81), if the desire is for all 0 o 
to be accelerated quickly. 

Let us once more consider a similar analysis of different G 0 but with the lowest an initial 
electron energy that will satisfy the resonance condition. We have y 0 = 7min and a res = 0 
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Figure 5.7: Simulation results for on-resonance but different initial 0 o for y () = 119, B„ = 45. The 
top graph shows energy versus distance traveled for the first 15mm. The middle graph 
shows /3 ± versus distance. The bottom graph shows the evolution of 6. The black curves 
shows d 0 = 71. The green curves show 6,, = n/2. Red curves show 0 o = 0. 


rad. These results are shown in Fig. 5.8. Minimum energy injection results in /3j_ = 0. To 
prove this, substitute the analytic expression for minimum energy required for resonance, 
Eq. (3.81), into the expression for /3j_ in terms of 7 and e, Eq. (3.29). With = 0, the 
vector j3 ± has no length. Therefore, 0 is undefined. In fact, all electrons starting at the 
minimum energy, will have the same evolution for energy and /3j_ as shown in Fig. 5.8. 
The laser accelerator for energies close to the minimum required for resonance will accept 
all values of 6 0 and resulting in energy gain. 


58 










150 



4 X1Q ^ 



Figure 5.8: Simulation results for on-resonance with y 0 = y mm ~ 118.99927, B„ = 45. The top graph 
shows energy versus distance traveled for the first 15mm. The bottom graph shows /i 
versus distance. 


Let us now try to understand why this surprising result makes physical sense. Recall that 
this entire discussion has been on resonance or when Aft) = 0. Recall Eq. (3.53) where 
tj = y(Aft) + ft>/,j3_sin 0). Since we know that is zero at some point in the evolution 
for the 0 o = 0 case, that pins the value of r/ to zero for all time since it is a constant of the 
motion. But later in the evolution, /3j_ and y grow. The only way to ensure that r/ stays zero 
is for sin 0 to be zero or when 0 — K. The next line of reasoning eliminates the possibility 
that 9 evolves to 0. 

When /3j_ is zero, then /3 V and /3 V are also zero. Going back to the original equations of 
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motion, Eqs. (3.19) through (3.22), we now have 


d -l = 0 

dt 

d 

eE n 

cos £ 

= - 

■ — ( 1 -/ 3 .) 

me 

// 

eE n - , 

sin£ 

= 

— 1 ~Pz) 

me 

4(rft) 




(5.4) 

(5.5) 

(5.6) 

(5.7) 


So the terms that deal with the magnetic field are completely turned off. They contribute 
nothing when the electrons move parallel to the field lines. Rewriting Eqs. (5.5) and (5.6) 
with the components of the electric field factored out we have 



= ( E x ) 

ymc v 

(5.8) 

>.v) 

- ^ ( Ey) 

ymc 

(5.9) 


The right-hand side of Eqs. (5.8) and (5.9) are proportional to the force that the electron 
feels when it is moving in a direction that is opposite to the electric field or when 9 — K. 

In this section we have considered the performance of the laser accelerator for only three 
9 a cases. But we are encouraged that even with the most unfavorable initial 9 0 , energy gain 
is possible if the initial energies are near the minimum required for resonance. 


5.3 Off-resonance Performance 

As already discussed, the electron bunch in an actual experiment is much longer than the 
optical wavelength and therefore spans over the entire 9 0 range from 0 to 2k. Additionally, 
electrons will not be perfectly collimated and will consequently have a finite spread in 
injection angles a. This latter effect implies that the resonance condition, Aft) = 0, will not 
be satisfied for all electrons in the bunch since for a given initial energy, magnetic field, 
and laser frequency, A ft) is a function of the injection angle a. Substituting Eqs. (3.24) and 
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(3.41) in (3.43) we have 


A£0 = W6 ~i 1 ~ Vo \frf-~ lcosa ) ■ (5 - 10) 

The right-hand side of Eq. (5.10) equals zero when a — a res (which is calculated in Eq. 
(3.79)). We want to explore the behavior of the laser accelerator for off-resonance cases for 
a spread of a near a res . 

Let us begin with the case y c = 130 with an external field B u = 45 T. The results are shown 
in Fig. 5.9 where we have conducted many simulations of the laser accelerator over a meter 
long interaction and plotted the final energy values Yf as a function of a and 9 a . Notice 
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Figure 5.9: Surface graph of final electron energy Yf versus initial 0„ and a for y„ = 130 and B (l 
45. There is a peak energy gain near a^s ~ 2.4 mrad and 6 a = n as expected. 


that there is a peak energy gain near a res ~ 2.4 mrad and 0 o — K as expected. 

We can get a better perspective by considering the contour plot shown in Fig. 5.10. The 
results indicate that electrons with initial conditions inside the blue region achieve 99 per- 
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Maximum y = 2056 



o (mrad) 

Figure 5.10: Contour graph of electron energy versus 9 a and injection angles a for y 0 = 130 and B„ 
= 45. The peak occurs near 0,-es ~ 2.4 mrad. Electrons with 0,, and a inside the blue 
region achieve 99 percent of maximum energy gain. Electrons with 6 a and a inside 
the green region achieve 99.1 percent of maximum energy gain. Electrons with 6 a and 
a inside the maroon region achieve 99.9 percent of maximum energy gain. 


cent of maximum energy gain. Electrons with initial conditions inside the green region 
achieve 99.1 percent of maximum energy gain. Electrons with initial conditions inside the 
red region achieve 99.9 percent of maximum energy gain. These results give a sense of 
the energy spread of the laser accelerator output beam, depending on the alignment of the 
input beam from the electron injector and the spread in injection angles. 

Let us consider a similar plot but now for the case of y () — 120 with an external field of 
45 T shown in Fig. 5.11. These results look similar to the previous contour plot but the 
parameter space area resulting in 99 percent energy gain has increased dramatically to 
include all initial 6 0 values. From the discussion on phase from § 5.2, this makes sense. 
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Maximum y = 2044 



o (mrad) 

Figure 5.11: Contour graph of electron energy versus 6 () and injection angles a for y 0 = 120 and B„ 
= 45. The peak occurs near Ores ~ 1.8 mrad. Electrons with 0„ and a inside the blue 
region achieve 99 percent of maximum energy gain. Electrons with 9 a and a inside 
the green region achieve 99.1 percent of maximum energy gain. Electrons with 9„ and 
a inside the maroon region achieve 99.9 percent of maximum energy gain. 


Electrons with an unfavorable 9 0 , lose /3_ since J3j_ begins with a smaller value. After /3_ 
goes to zero, The value of 9 is no longer meaningful and can change to % for the reasons 
discussed previously. Now these electrons now are able to achieve energy gain. 

Let us consider a similar plot but now for the case of y 0 = 119 with an external field of 45 
T shown in Fig. 5.12. The parameter space area resulting in 99.1 percent energy gain has 
increased dramatically to include all initial 9 0 . Notice that for Figs. 5.10, 5.11, and 5.12, 
we gained greater acceptance range in 9 0 at the expense of the maximum energy j/ since 
it went down from jf — 2056 to jf = 2042. 
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Maximum y = 2042 



o (mrad) 

Figure 5.12: Contour graph of electron energy versus 9 a and injection angles a for y 0 = 120 and B„ 
= 45. The peak occurs near Ores ~ 20 /.(rad. Electrons with 8 () and a inside the blue 
region achieve 99 percent of maximum energy gain. Electrons with 6 a and a inside 
the green region achieve 99.1 percent of maximum energy gain. Electrons with 6 a and 
a inside the maroon region achieve 99.9 percent of maximum energy gain. 


Encouraged by the previous contour plots and the results of § 5.2, it seems that starting 
with a smaller /i L ensures more electrons receive similar accelerations and energy gains. 
So let us just start with a beam entering on-axis where a — (5^ — 0. Now Aft) is only a 
function of two independent variables, y 0 and B u . We have for on-axis injection 

A®oa =--0)r(l-—Jy$-l). (5.11) 

my 0 V Yo v ) 

By varying y () and B„ we expect the maximum energy to be in a band that follows the 
A ft) = 0 contour seen in Fig. 5.13. 
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Figure 5.13: Simulation results showing an intensity plot of final electron energies with initial con¬ 
ditions that explore a range of external magnetic fields from B„ = 30 T to 60 T and 
y 0 = 80 to 130. The dotted line indicates the Acoqa = 0 contour. 


In Fig. 5.13, we have an intensity plot that depicts the final electron energies of electrons 
with varying initial y 0 and external magnetic field values but with on-axis injection where 
a = /3j_ = 0. Therefore the angle between j3 ± and E r has no meaning. The optical force 
on the electron from the electric field acts on the electron such that 0 quickly evolves to 
a value of n as discussed at the end of § 5.2. The electrons then proceed to gain energy. 
The resonance A (o = 0 level curve is overlaid on the intensity plot and shown as a black- 
dotted curve. Notice how the electrons with good energy gain follows that contour curve as 
expected. There exists a width about the contour line A a> = 0 for whereby the laser accel¬ 
erator performs well with jj « 2500. This corresponds to electron energy gain of 1.2 GeV 
in one meter. Once the right combination of parameters for resonance are identified, the 
laser accelerator is quite forgiving in its acceptance margins in initial energy and external 
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magnetic field strength for on-axis injection. 

We can now make an estimate on normalized emittance values that the laser accelerator 
would accept from an electron injector. Normalized emittance is given by 

£ n = Y<r><a> (5.12) 

Let us consider an initial electron energy of y = 130, < r > =0.5 mm since the optical waist 
is about 1 mm, and < a > = 0.5 mrad. The spread in injection angles a is a conservative 
estimate based on the simulations presented. We then have £ n = 32.5 mm-mrad, which is 
quite large. Electron injectors with good emittance values are about 10 mm-mrad. 


5.4 Magnetic Field Effects 

The external magnetic field source for the laser accelerator might be constructed of resistive 
or superconducting solenoids such as those at NHMFL. The magnetic field on the axis for 
a solenoid can be written as [25]: 


B- 5 m (cos0i + cos 62 ) (5.13) 

where 0\ and 02 are defined in Fig. (5.15). 

For example, if we assume the maximum field 5™ ax = 45 T at z w — 0.5 m, then by Eq. 
(5.13) the field profile is given by Fig. (5.14) assuming a bore diameter of 32 mm. 

By including Eq. (5.13) in the simulation we have the following results shown in Fig. 5.16. 
Notice the results indicate a final energy of 1.15 GeV versus 1.2 GeV for a the static case. 
The top left graph plots electron energy gained versus distance traveled. The bottom left 
graph plots /3j_ versus distance traveled. The bottom right graph plots 1 — j3 z versus distance 
traveled. The only major difference with Fig. 5.2 is that the electron’s radial position flared 
out in the last few millimeters. Comparing this with Fig. 5.15, the magnetic field strength 
drops to nearly half its maximum value. Therefore, the electron spirals out more. 
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Figure 5.14: Magnetic field variation along z with maximum field at z w = 0.5 m [25] 



Figure 5.15: Diagram defines 0 1 and 6F for calculating the magnetic field along the longitudinal 
axis of a solenoid represented by the box. 


5.5 Diffraction Effects 

The numerical solutions presented in the previous sections have assumed plane-wave de¬ 
scription for the laser source. In real lasers, the optical beam wavefronts are spherical with 
the fundamental mode taking on a Gaussian irradiance pattern transverse to the optical axis 
. fields for the electric field of the laser where the wavefronts ar. Real lasers have a Gaussian 
approximation to the laser However, let us consider the effects of Gaussian mode. Begin- 
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Figure 5.16: Real magnetic field effects included in simulation for y a = 119. The top left graph 
plots electron energy gained versus distance traveled. The top right graph plots the 
radial position of the electron from optical axis. The bottom left graph plots f5± versus 
distance traveled. The bottom right graph plots 1 — J3- versus distance traveled. 


ning with a circularly-polarized Gaussian laser field in the fundamental mode, we have the 
following electric and magnetic field components: 


E r 


^ e ~r 2 /w 2 (z) e WH) 

w(z) 



! r = ^^ e ~> 2 /w 2 (z) e i^H) 

1 c w{z) 1 


(5.14) 
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with 


, _ k, r 2 

^ G ~ 2 R 

"(V) 

(5.15) 

w = w 0 J 

\ . (z-z w ) 2 

z 2 

(5.16) 

nwl 



y 0 

Lr ~ A 


(5.17) 

R = - i —[(z-z w ) 2 + Z;] 

Z — Zw 

(5.18) 


where £ =k r z — CO r t, (pc is the Gaussian phase contribution, w is the beam radius radius, 
w a is the beam waist, Z, is the Rayleigh length, R is the phase front radius of curvature, and 
z w is the position of the laser beam waist. 

Substituting Eqs. (5.14) into (3.4) then substituting the real part into Eq. (3.3) yields the 
following: 

dj_ = |ftcos(fc + f) _ j3vSin(fc+ f)| (5 . 19 ) 

dt me w 

d , 0 . eE a Wo _ r 2/ w 2 n \ _ 

tMPx) = - e ' COS0G + C 1-A- Av (5-20) 

dt me w m 

4(jft,)= — sin(fc + {)(1 - ft) + —ft (5.21) 

at me w m 

4(rft) = [ftcosffc + 0 - ft sin(fc + 0] (5.22) 

dt me w 

Eqs. (5.19) through (5.22) are similar to Eqs. (3.19) through (3.22) in the previous section 
but they are now modified to include Gaussian beam effects in the amplitude and phase 
terms. To first order, let us assume the following: the electron’s trajectory is well inside 
the optical mode and the Rayleigh length, Z r , is two meters long (twice the interaction 
region distance). With these assumptions, the amplitude factor, w 0 e~ r ' w /w, does not 
significantly change since the electron radius is well inside the optical mode. The Guoy 
phase shift, (pc, is 

(j) G = - tan -1 • (5-23) 
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Figure 5.17: Guoy phase shift over the interaction region with z w = 0.5. 


and plotted in Fig. (5.17). As we expect, to first order, the d(j) G /dz ~ — 1 jZ r about the 
waist. 

Once again we can define the constant of the motion, £ — y( 1 — /3-), based on Eqs. (5.19) 
and (5.22) by noting that the right-hand sides are identical. Now we must consider how (pc 
affects the other constant of the motion, 77 . 

Note the following: 

-^cosOg+C) = -sin(0G + C)^ + ®r(l-/3 z )sin(^G + C) (5-24) 

dt dt 

^sin(0 G + O = cos(<^G + C)^-®r(l-/3 z )cos(<^G + C) (5-25) 

sin(0G + O^(yAr) = J t (rA^in(0c + 0) - Yl^xcos((j) G + + YPxM 1 - &)cos((^G + C) 

(5.26) 

cos(0G + C)^(rAy) = J t (rAyCOs(0G+C)) + yAySin(0G + O^-rAy®r(l -Az)sin(0G + O 

(5.27) 
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Solving for /3 Y and /3 y in Eqs. (5.20) and (5.21) and substituting in Eq. (5.19) we have 


dy 

dt 


m 1 a sdy 

yen, 1-/3- —- 

e B„ dt 


r 


d<j)c m dy E„ d 


dt eB u dt 


cB u dt 


(rAv'COS(0G + O)- 


Eo d 
cB u dt 


(y/3 x sin(0c + C))- 

(5.28) 


To first order, dfyojdt = (d^c/dz) ( dz/dt ) ~ — c/Z r . Substituting this into Eq. (5.28) and 
extracting an overall time derivative we arrive at 


eBu 

y 0 m 


®r(l ~Pzo) H-- (/^a- sin(0c + C) + A’ c os(0c+ 0) 

y 0 mc 


c y 
2Z, y, 


= constant 
(5.29) 


and as before Acu = eB u /y 0 m — co, ( 1 — j8 z0 ), y 0 = y(0), e = y 0 ( 1 — /3- 0 ), and 0 )e — eE () /ymc. 
Comparing Eq. (5.29) to Eq. (3.44) we see that the Gaussian beam description has added 
0c in the trigonometric functions and a new term, —cy/2Z r y 0 . The magnitude of this first 
order correction is on the order of 10 9 s^ 1 . This is approximately two-orders of magnitude 
smaller than other terms. The conclusion is that, to first order, a Gaussian beam will have 
similar solutions to the plane-wave solutions calculated in the previous sections. 
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CHAPTER 6: 
Conclusions 


The relativistic Lorentz force law for the laser accelerator has been solved numerically 
using Euler-Cromer integration methods for the first time. Simulation results indicate that 
electron energy can grow from about 60 MeV to over 1.2 GeV in one meter using today’s 
high powered lasers and magnetic fields. The electron’s perpendicular velocity is a key 
parameter in achieving good energy gain. For the first time, the complete evolution of 
the electron’s perpendicular velocity over the interaction region is shown. The resonance 
condition was reexamined and a minimum energy requirement was determined for the first 
time. Laser accelerator performance was examined when at A co / 0 and was determined 
to be more tolerant of off-resonance injection than previously thought. 

The initial electron phase, or the angle between /3 ^ and E r , was examined and surprisingly 
found to be tolerant of unfavorable initial phase orientation. If operating near the minimum 
energy and at resonance, it was found that there is significant energy gain for all initial 
electron phases. Electrons that were once thought to not gain energy due to unfavorable 
initial phase orientations may still contribute to a high energy beam if the initial energy is 
near the minimum energy. Therefore, a large parameter space in phase and injection angles 
exists for the laser accelerator. 

Diffraction, magnetic field solenoid edge, and Guoy phase evolution effects were examined 
and differences from ideal plane wave results are provided. Magnetic field edge effects 
were examined and determined to not adversely affect electron energy gain. 

A survey of existing high powered laser sources, magnetic field systems, and injectors was 
conducted to give a sense of what a university sized laser accelerator would have to be. 
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APPENDIX: Simulation Codes 


There were many individual programs that were created in C or Matlab for a specific task 
such as calculating resonance angle for a given magnetic field, laser wavelength, and initial 
electron injection energy. In this appendix a sample of the laser accelerator engine code 
is included for reference. The numerical integration was programmed in C with graphics 
exported for use in Matlab. 

A.l Laser Accelerator Simulation Core 

// 

// main.c 
// LAsim7 
// 

// Created by Ricardo Vigil on 3/4/14. 

// Copyright (c) 2014 Ricardo Vigil. All rights reserved. 

// 

#include <stdio.h> 

#include <stdlib.h> 

#include <math.h> 

#include <time.h> 

#define pi 3.1415926535897932 

//double pi = 3.141592653589793238; 

double c = 2.99792458E8; 
double uo = 12.566370614E-7; 
double eo = 8.854187817E-12; 
double q = 1.602176565E-19; 
charge (C) 

double m = 9.10938291E-31; 


// speed of light (m/s) 

// Permeability of vacuum (N A A -2) 
// permittivy constant (F/m) 

// absolute value of electron 

// mass of electron (kg) 
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double mc2_ev = 0.510998928E6; 

double Zr = 2.0; 

double lambda_r = IE-6; 

laser 

double k; 

double P = 1E15; 

is 10 A 15 

double Ipk = 1.0E25; 

Vulcan Laser 

double delta_z = IE-5; 

size (m) 

double L = 1.0; 

double wo; 

m for 

double xo = 0.0; 
double yo = 0.0; 
double dt = 7E-14; 
double tiempo; 
double Bu = 45.0; 
double width; 
double wr; 
double epsilon; 
double dw; 

double gamma_min, gamma_o, ggamma; 
double alpha; 

double A, Eo, beta_o, we, bzo, bpo 
g_star; 

double zeta, bp, bz, bx, by, bz, x 
g_check_bxby; 

double ax, ay, az, dx, dy, dz, t, 
double bxold, byold, bzold, LP; 
double phi_e, phi_l, bxl, byl, Ex, 


// elctron rest mass (eV) 

// Rayleigh length (m) 

// laser wavelength (m) Vulcan 

// wave number (m A -l) 

// laser peak power (W) source 

// peak intensity W m A -2 at 

// propagation length step 

// length of interest (m) 

// laser beam waist radius in 

// initial electron position 
// initial electron position 
// time step 

// time for light to travel 
// External B field 
// pipe width (m) 


wco, g_m, bpmax, rot_angle, 
Yr z, g_check_bp, 

, a,b,pi,p2,bxo, byo; 

Ey, brk, Bu, r_approx. 
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alpha_ 

double 

dalpha 

alpha_ 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 

double 


o, phi, eta; 

gamma_top, gamma_bot, dgamma, alpha_top, alpha_bot, 

, epsilon_res, 
res, kappa; 

pO, px, py, pz, ggammaold, Eng, eta, bp_bxby, bp_epsilon; 
pOold, pxold, pyold, pzold; 
ctheta, stheta, stheta2, brk, brk2; 
dl[1000001] ; 

delta_g, g, intG, epsilon_f; 
frac_a, frac_b, frac_ep, frac_eta; 
ao, bo, epsilon_o, eta_o, dwo; 

Ex, Ey; 

Eangle, bpangle; 
etap, etam; 
theta; 

etaprime, temp, wtf; 
bl, b2; 

ggammal, zl, test; 
aa, bb, cc, zz; 
zo; 

frac_dw, dw_anal, dw_te, alpha2, gp, gm; 
denom, denom_ul, denom_ul_explicit; 
factorl, factor2, factorl2, denomFactor; 
factor3, factor3eta; 
factor4; 

expandedFactor3; 
expandedDenom; 
expandedDenom_dwwe; 
a a_u, bb_u, c c_u; 

d, f; 

dwwe_zero, dwwe_l, dwwe_2, LHS, RHS; 
delta, rho; 
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double deltaomega; 
double z_anal; 


int i,j; 
int nalpha=150; 
int ngamma=lll; 
int nt; 

int count=1000001; 


int flag = 4; 


// 

case 

1. 

resonance & 

consts 

// 

case 

2 . 

resonance & 

no 

consts 

// 

case 

3. 

no resonance 

& 

consts 

// 

case 

4 . 

no resonace 

& 

no consts 

// 

case 

5 . 

explict use 

of 

cos(theta 


int main(void) 

{ 

clock_t begin, end; 
double time_spent; 
begin = clock(); 

FILE *out10, *outll; 
FILE *out30, *out31; 


78 



// - 

// Initialize arrays 

// - 

for(i=0;i<1000001;i++) dl[i] = 0.0; 


// - 

// Calculate Global constants 

II - 

nt = round(L/(c*dt)) + 1; 

// number of time steps 
width = 32.0 / 1000; 

// width of pipe. 

wo = sqrt(Zr*lambda_r/pi); 

// mode radius at waist 
A=pi*wo*wo; 

// area of laser beam (m A 2) at waist 
Eo = sqrt( (P)/(c*eo*A)); 

// Electric Field (V/m) 

wr = (2*pi*c/lambda_r); 

// laser field angular freq 
k = (2*pi)/lambda_r; 

// laser wave number 
deltaomega = 0.0; 

printf("wr\t\t\t = %.16e\n",wr); 

printf("deltaomega\t = %.16e\n",deltaomega); 

wr = wr + deltaomega/2.0; 

printf("wr (new)\t = %.16e\n",wr); 


// 
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// Initial Conditions (Global) 

// - 

printf("Eo = %.16e\n",Eo); 

x = 0.0; 
y = 0.0; 
z = 0.0; 
z_anal = 0.0; 
r = x*x + y*y; 
r_approx = 0.0; 
t = 0.0; 

LP = 0.0; 
rot_angle = 0.0; 
zeta = k*z - wr*t; 

switch (flag){ 

case 4 : 

{ 

printf("case 4.\n\n"); 
printf ( "wo\t\t = %.16e\n",wo); 

FILE *out40, *out41, *out42, *out43, *out44; 

Bu = 45.0; 

out40 = fopen("/Users/vigilr/Dropbox/FEL/ 
LAsim7/lasim7_case4_Io.txt","w"); 

out41 = fopen ( "/Users/vigilr/Dropbox/FEL/ 
LAsim7/Iasim7_case4.out","w"); 

out42 = fopen("/Users/vigilr/Dropbox/FEL/ 
LAsim7/Iasim7_case4_int.out","w"); 

out43 = fopen("/Users/vigilr/Dropbox/FEL/ 
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LAsim7/Iasim7_case4_test.out","w"); 

out44 = fopen("/Users/vigilr/Dropbox/FEL/ 
LAsim7/Iasim7_case4_fracDW.out","w"); 


// - 

// Initial Conditions (specific to Case 4) 
//No constraints on gamma_o, alpha, 
and phi for a fixed Bu. 

// - 

ggamma = 150.0; 
ggammaold = ggamma; 
gamma_o = ggamma; 

phi = pi; //3.1292756331662902e-01; 
alpha_res = acos((gamma_o - ((q*Bu)/(wr*m)) )/ 

(sqrt(gamma_o*gamma_o - 1)) ); 

alpha = alpha_res; 


bz = 1.0/ggamma*sqrt(ggamma*ggamma - 1.0)*cos (alpha); 
printf("bzo\t\t= %.16e\n",bz); 

bp = 1.0/ggamma*sqrt(ggamma*ggamma - 1.0)*sin(alpha); 
printf("bpo\t\t= %.16e\n",bp); 

bx = bp*cos (phi); // eq 4 

by = bp*sin(phi); // eq 5 

beta_o = 1.0/ggamma*sqrt(ggamma*ggamma - 1.0); 
if (fabs (bp) < IE-15) bp=0.0; 
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if (fabs(bx) < IE-15) bx=0.0; 
if (fabs (by) < IE-15) by=0.0; 
zo = atan (bx/by)/k; 

epsilon = ggamma*(1-bz); // eg 6 


bp_bxby = sqrt(bx*bx + by*by); 
bp_epsilon = (1/ggamma)* 
sqrt(2*ggamma*epsilon - epsilon*epsilon - 1); 


bpo = 

bp_epsilon; 




bxold 

= bx; 





byold 

= by; 





bzold 

= bz; 





pO = 

ggamma*m*c; 

// 

eq 

10 

px = 

ggamma*m*c*bx; 

// 

eq 

7 

py = 

ggamma*m*c*by; 

// 

eq 

8 

pz = 

ggamma*m*c*bz; 

// 

eq 

9 

pOold 

ii 

PS 

o 





pxold 

= px; 





pyold 

ll 

P 

P 





pzold 

= pz; 





g_check__bxby = 

1/sqrt(1 

- bx*bx - 

b 

i 

b 

-X 

ax = 

c*bx/dt; 





ay = 

c*by/dt; 





az = 

c*bz/dt; 





dx = 

c*by*az - 

ay*c*bz; 




dy = 

ax*c*bz - 

c*bx*az; 




dz = 

c*bx*ay - 

ax*c*by; 




zeta 

o 

o 

II 

*z - wr*t; 




Ex = 

cos(zeta) , 

r 




Ey = 

-1.0*sin (: 

zeta); 




fprintf(out43,' 

"%.16e\n", 

(pxold*cos(z 

eta 


- bz*bz); 


- pyold * 


// 


sin(zeta) 
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//- 

— 

— 

— 

— 

— 

— 

— 

— 



// 


initial angles 




//- 

— 

— 

— 

— 

— 

— 

— 

— 



if 

(Ex 

> 

0 

&& 

Ey 

> 

0 ) 

Eangle = 

atan(Ey/Ex); 


if 

(Ex 

< 

0 

&& 

Ey 

> 

0 ) 

Eangle = 

atan(Ey/Ex) 

+ pi; 

if 

(Ex 

< 

0 

&& 

Ey 

< 

0 ) 

Eangle = 

atan(Ey/Ex) 

+ pi; 

if 

(Ex 

> 

0 

&& 

Ey 

< 

0 ) 

Eangle = 

atan(Ey/Ex) 

+ 2.0 *pi; 

if 

(bx 

> 

0 

&& 

by 

> 

0 ) 

bpangle 

= atan(by/bx) 

f 

if 

(bx 

< 

0 

&& 

by 

> 

0 ) 

bpangle 

= atan(by/bx) 

+ pi; 

if 

(bx 

< 

0 

&& 

by 

< 

0 ) 

bpangle 

= atan(by/bx) 

+ pi; 

if 

(bx 

> 

0 

&& 

by 

< 

0 ) 

bpangle 

= atan(by/bx) 

+ 2.0 *pi ; 


theta = Eangle - bpangle; 

if (theta < 0 ) theta = theta + 2.0*pi; 

if (theta > 0 && theta > pi) theta = 2.0*pi - theta; 

if (phi>0 && phi<pi) theta = phi; 

// - 

// Calculate derived paramaters 
//- 

test = bx*cos (k*z-wr*t+k*zo) - by*sin(k*z-wr*t+k*zo); 

we = (q*Eo)/(gamma_o*m*c); 

wco = (q*Bu)/(m*ggamma); 

g_star = (epsilon*epsilon+l)/epsilon; 

gamma_min = (epsilon*epsilon+l)/(2*epsilon); 

bpmax = (epsilon/sqrt(epsilon*epsilon+l)) 

dw = wco + 

(wr/ggamma)*sqrt(ggamma*ggamma - 1)*cos(alpha) - wr; // eq 23 
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dw_anal = (q*Bu)/(ggamma*m) + 

wr*bz - wr; 

a = ggamma*bx - ((q*Eo)/(m*c*wr))*sin(zeta) 

+ ((q*Bu)/(m*c))*y; 

b = ggamma*by - ((q*Eo)/(m*c*wr))*cos(zeta) 

- ((q*Bu)/(m*c))*x; 

temp = (bx*sin(zeta) + by*cos(zeta)); 
eta = ggamma*dw + ggamma*we*(bx*sin(zeta) 

+ by*cos (zeta)); // eq 24 

etap = ggamma*dw + 
ggamma*we*bp_epsilon*sin(theta) ; 

etam = ggamma*dw - 
ggamma*we*bp_epsilon*sin(theta) ; 

wtf = asin((wr - wco - wr*beta_o*cos(alpha)) 

/(we*beta_o*sin(alpha))) ; 

stheta = (eta - ggamma*dw)/ 

(we*sqrt(2*ggamma*epsiIon - epsilon*epsiIon - 1)); 

stheta2 = (bx*sin (zeta) + by*cos(zeta))/bp_bxby; 
ctheta = (bx*cos (zeta) - by*sin(zeta))/bp_epsilon; 
brk = bx*cos(zeta) - by*sin(zeta) ; //bp * cos (theta); 
brk2 = bx*sin(zeta) + by*cos(zeta); // bp * sin(theta); 

frac_a = 0.0; 
frac_b = 0.0; 
frac_ep = 0.0; 
frac_eta = 0.0; 

cLO = 3l} 

bo = b; 

epsilon_o = epsilon; 
eta_o = eta; 

//- 

// Inputs 
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// - 

fprintf(out40,"case 4. no resonance & consts\n"); 
fprintf (out40, "\t-Inputs-\n") ; 


// - 

// Equations of Motion 

II - 


fprintf(out41,"%.16e\t%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\t%.16e\t%16\n", 
t,x,y,z,r,r_approx,ggamma,bx,by,bz, 
bp_bxby,bp_epsilon,g_check_bxby,LP,zeta, 
stheta,ctheta,brk,epsilon,a,b,eta, 
frac_a,frac_b,frac_ep,frac_eta, 
brk2,stheta2,theta,etap,etam, 
etaprime,dw_anal,z_anal); 

for(i=0;i<nt;i++){ 

//— 
pO = pO + 

((-q*Eo/(ggammaold* m*c))*(pxold*cos(zeta) - pyold * sin(zeta) ) 

ggamma = pO/(m*c); 

bz = (ggamma - epsilon)/ggamma; 

z = z + bz*c*dt; 


) *dt; 
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//- 

px = pxold + ((-q*Eo)*(epsilon/ggammaold)*cos(zeta) 

- (q*Bu*pyold)/(ggammaold*m))*dt; // eq 15 
bx = px/(ggamma*m*c); 
x = x + bx*c*dt; 

//- 

PY = pyold + ((q*Eo)*(epsilon/ggammaold)*sin(zeta) 

+ (q*Bu*pxold)/(ggammaold*m))*dt; // eq 18 
by = py/(ggamma*m*c); 
y = y + by*c*dt; 

//- 

pxold = px; 
pyold = py; 
ggammaold = ggamma; 

//- 

t = t + dt; 

zeta = k*z - wr*t; 

//- 

bp_bxby = sqrt(bx*bx + by*by); 

bp_epsilon = (1/ggamma)*sqrt(2*ggamma*epsilon 
- epsilon*epsilon-l); 

if ((2*ggamma*epsilon - epsilon*epsilon-l) < 0 ) { 

bp_epsilon = 0.0; 

} 

//- 

stheta = (eta_o - ggamma*dw)/ 

(we*sqrt(2*ggamma*epsilon - epsilon*epsilon - 1));; 

stheta2 = (bx*sin(zeta) + by*cos(zeta))/bp_bxby; 
ctheta = (bx*cos(zeta) - by*sin(zeta))/bp_epsilon; 
brk = bx*cos(zeta) - by*sin (zeta); // 
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brk2 = bx*sin(zeta) + by*cos(zeta); 
epsilon_f = ggamma*(1-bz); 


//- 

— 

— 

- 

— 

— 


— 


// 


Angles 





//- 

— 

— 

- 

— 

— 


— 


Ex 

= 1 

SOS 

( 

zeta); 




Ey 

= 

sin ( 

zeta); 




if 

(Ex 

> 

0 

&& 

Ey 

> 0 

) Eangle = 

atan(Ey/Ex) ; 

if 

(Ex 

< 

0 

&& 

Ey 

> 0 

) Eangle = 

atan(Ey/Ex) 

if 

(Ex 

< 

0 

&& 

Ey 

< 0 

) Eangle = 

atan(Ey/Ex) 

if 

(Ex 

> 

0 

&& 

Ey 

< 0 

) Eangle = 

atan(Ey/Ex) 

if 

(bx 

> 

0 

&& 

by 

> 0 

) bpangle 

= atan(by/bx) 

if 

(bx 

< 

0 

&& 

by 

> 0 

) bpangle 

= atan(by/bx) 

if 

(bx 

< 

0 

&& 

by 

< 0 

) bpangle 

= atan(by/bx) 

if 

(bx 

> 

0 

&& 

by 

< 0 

) bpangle 

= atan(by/bx) 


theta = Eangle - bpangle; 

if (theta < 0 ) theta = theta + 2.0*pi; 

if (theta > 0 && theta > pi) theta = 2.0*pi 

bl = bp_epsilon*cos(theta); 
b2 = bp_epsilon*sin(theta); 

//— 

a = ggamma*bx - ((q*Eo)/(m*c*wr))*sin(zeta) 

+ ((q*Bu)/(m*c))*y; 


pi; 

pi; 

2.0 *pi; 

+ pi; 

+ pi; 

+ 2.0 *pi 


theta; 
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b = ggamma*by - ((q*Eo)/(m*c*wr))*cos (zeta) 

- ((q*Bu)/(m*c))*x; 

eta = ggamma*(dw + ((q*Eo)/(gamma_o*m*c))* 

(bx*sin (zeta) + by*cos(zeta))) ; 

etap = ggamma*(dw + ((q*Eo)/(gamma_o*m*c))* 
bp_epsilon*sin(theta) ) ; 

etam = ggamma*(dw - ((q*Eo)/(gamma_o*m*c))* 
bp_epsilon*sin(theta)); 

etaprime = ggamma + ggamma*(Eo/c*Bu)* 

(bx*sin (zeta) + by*cos(zeta)) - (m*epsilon*wr*ggamma)/(q*Bu) 

frac_a = (a - ao)/ao; 

frac_b = (b - bo)/bo; 

frac_eta = (eta - eta_o)/eta_o; 

frac_ep = (epsilon_f - epsilon)/epsilon; 

// — 

r = sqrt(x*x + y*y); 

r_approx = (m*ggamma*sqrt(bx*bx+by*by)*c)/(q*Bu); 

//- 

ax = c* (bx - bxold)/dt; 
ay = c*(by - byold)/dt; 
az = c* (bz - bzold)/dt; 
bxold = bx; 
byold = by; 
bzold = bz; 

dx = c*by*az - ay*c*bz; 

dy = ax*c*bz - c*bx*az; 

dz = c*bx*ay - ax*c*by; 

LP = ((uo*q*q*ggamma*ggamma*ggamma*ggamma* 

ggamma*ggamma)/(6*pi*c) )* 

(ax*ax + ay*ay + az*az - 

(1/ (c*c))*(dx*dx + dy*dy + dz*dz)); 


88 



//- 

rot_angle = rot_angle + (q*Bu)/(m*ggamma)*dt; 

//- 

g_check_bxby = 1/sqrt(1 - bx*bx - by*by - bz*bz); 
//- 

dw_anal = (q*Bu)/(ggamma*m) - wr*(1-bz); 

// - 

z_anal = (c/gamma_o)* 

(1/ (3*we*we*we*we*epsilon*epsilon) ) * 
sqrt(2*we*we*epsilon*ggamma - 
we*we*epsilon*epsilon - we*we - eta*eta)* 
(we*we*epsilon*ggamma - 
we*we*epsilon*epsilon + we*we + 
eta*eta - 3*we*we*epsilon*epsilon) ; 


// - 

// Matlab output 
//- 


fprintf(out41,"%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\t%.16e\t%.16e\t 
%.16e\t%.16e\n", 
t,x,y,z,r,r_approx,ggamma, 
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bx,by,b z,bp_bxby, 

bp_epsilon,g_check_bxby,LP, 

zeta,stheta,ctheta,brk,epsilon_f, 

a,b,eta,frac_a,frac_b,frac_ep,frac_eta, 

brk2,stheta2,theta,etap,etam,etaprime,dw_anal,z_anal); 
} 

// - 

// Outputs 

// - 

epsilon_f = ggamma*(1-bz); 
a = ggamma*bx - 

( (q*Eo)/(m*c*wr) )*sin (zeta) + ((q*Bu)/(m*c))*y; 

b = ggamma*by - 

( (q*Eo)/(m*c*wr) )*cos (zeta) - ((q*Bu)/(m*c))*x; 

eta = ggamma*(dw + 

((q*Eo)/(ggamma*m*c))*(bx*sin(zeta) + by*cos(zeta))); 


fprintf(out40,"\n\t-Calculated Values-\n"); 

fprintf(out40,"?_max\t= 

%.16e\n?pmax\t= %.16e\n",g_m,bpmax); 
fprintf (out40, 

"\n\t-Outputs-\n") ; 

fprintf(out40,"?_F\t= %.16e\n 
Eng_F\t= %.16e\n?z\t= %.16e\n 
?p\t= %.16e\nRev\t= %.16e\n 
r_F\t= %.16e\nr_apprx\t= %.16e\n 
L\t= %.16e\nTime\t= %.16e\n", 
ggamma,mc2_ev*ggamma/1E9,bz. 
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bp,rot_angle/(2*pi),r,r_approx,z,t); 

print f ( " ?_star\t= % . 16e\n?__min\t=% . 16e\n 
?_o\t\t= %.16e\n?_F\t\t= %.16e\n 
Eng_F\t= %.16e\n?\t\t= %.16e\n 
?\t\t= %.16e\nBu\t\t= %.16e\n 
?\t\t= %.16e\na\t\t= %.16e\n 
b\t\t= %.16e\n?\t\t= %.16e\n 
we\t\t= %.16e\nwc\t\t= %.16e\n 
wr\t\t= %.16e\n?w\t\t= %.16e\n 
Revs\t= %g\nnt\t\t= %d\nz\t\t= %.16e\n", 
g_star,gamma_min,gamma_o, 
ggamma,mc2_ev*ggamma/IE9,alpha, 
phi,Bu,epsilon,a,b,eta_o, 
we,wco,wr,dw,rot_angle/(2*pi),nt,z); 
fclose(out40); 
fclose (out41); 
fclose (out42); 
fclose(out43); 
fclose (out44); 
break; 

} 

default: 

printf("try again\n"); 

} // end switch 

// - 

end = clock (); 

time_spent = (double)(end - begin) / CLOCKS_PER_SEC; 
printf("\n\nElapsed: %.16g seconds\n", time_spent); 

return 0; 
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} 


A.2 Laser Accelerator Matlab Graphics Script 

clear 
close all 
clc 

load Iasim7_case4.out; 
load Iasim7_case4_int.out; 
load Iasim7_case4_lOm.out; 

c = 2.99792458E8; % 

uo = 12.566370614E-7; % 

eo = 8.854187817E-12; % 

q = 1.6 0217 65 65E-19; % 

m = 9.10938291E-31; % 

mc2_ev = 0.510998928E6; % 

Eo = 4.3401054911612872e+ll; 
we = (q*Eo)/(119.0*m*c); 


speed of light (m/s) 

Permeability of vacuum (N A A -2) 
permittivy constant (F/m) 
absolute value of electron charge (C) 
mass of electron (kg) 
elctron rest mass (eV) 

% 1/s 


%zl0 = Iasim7_case4_l0m(:,4); 
%ggammal0 = Iasim7_case4_l0m(:,7); 


t = Iasim7_case4(:,1); 
x = Iasim7_case4(:, 2) ; 
y = Iasim7_case4(:, 3) ; 
z = Iasim7_case4(:, 4) ; 
r = Iasim7_case4(:, 5) ; 
r_a= Iasim7_case4(:,6); 
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ggamma = Iasim7_case4(:,7); 
bx = Iasim7_case4(:,8); 
by = Iasim7_case4(:,9); 
bz = Iasim7_case4 (:,10); 
bp_bxby = Iasim7_case4(:,11); 
bp_epsilon = Iasim7_case4(:,12); 
g_check_bxby = Iasim7_case4(:,13); 
LP = Iasim7_case4(: , 14) ; 
zeta = Iasim7_case4(: , 15) ; 
stheta = Iasim7_case4(: , 16) ; 
ctheta = Iasim7_case4(: , 17) ; 
brk = Iasim7_case4(:,18); 
epsilon = Iasim7_case4(:,19); 
a = Iasim7_case4(:, 20) ; 
b = Iasim7_case4(:, 21) ; 
eta = Iasim7_case4(: , 22) ; 
frac_a = Iasim7_case4(: , 23) ; 
frac_b = Iasim7_case4(:,24); 
frac_ep = Iasim7_case4(:,25); 
frac_eta = Iasim7_case4(:,26); 
brk2 = Iasim7_case4(: , 27) ; 
stheta2 = Iasim7_case4(: , 28) ; 
theta = Iasim7_case4(: , 29) ; 
etap = Iasim7_case4(: , 30) ; 
etam = Iasim7_case4(: , 31) ; 
etaprime = Iasim7_case4(:,32) ; 
dw_anal = Iasim7_case4(: , 33) ; 
z_anal = Iasim7_case4(:,34); 

D = z_anal - (z_anal(2)-z_anal(1)); 
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zl = 

lasim7_ 

_case4 

_int ( : 

gl = 

lasim7_ 

_case4 

_int ( : 

ao = 

a 

(1) ; 



bo = 

b 

(i); 



f rac_ 

_a 

= (a 

- ao) 

/ao; 

frac_ 

_b 

= (b 

- bo) 

/bo; 


figure (1) 

subplot (3,2,1) 

set(gca,'FontSize',24) 

plot(z,ggamma*mc2_ev/IE9,' LineWidth' 

%hold on 

%plot(z,g_check_bxby*mc2_ev/IE9,' r', 

xlabel(' z','FontSize',24) 

ylabel('\gammamc A 2 (GeV)' ) 

xlim([0 max(z)]) 

grid 

legend('gamma','gamma check') 

subplot (3,2,2) 

%plot(z,r,'LineWidth',3) 
set(gca,'FontSize',24) 

%hold on 

plot(z,r_a,'r','LineWidth',2) 
xlabel ('z (m)') 
ylabel('r (m)') 
grid 

xlim([0 max(z)]) 


,3) 

' LineWidth',2) 
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legend(' r 


r approx') 


subplot(3,2,3) 

set (gca,'FontSize',24) 

plot(z,bp_epsilon,'LineWidth',3) 

%hold on 

%plot(z,bp_bxby,'r','LineWidth' , 2) 

%hold on 

%plot(z,brk,'g') 

xlabel ('z') 

ylabel('\beta\perp') 

xlim([0 max(z)]) 

grid 

subplot (3,2,4) 

set(gca,'FontSize',24) 

plot(z,1-bz,'LineWidth',3) 

xlabel ('z') 

ylabel('1 - \beta_z') 

grid 

xlim([0 max(z)]) 

subplot (3,2,5:6) 

plot(z,stheta,'LineWidth',3) 

hold on 

plot(z,ctheta,'r','LineWidth',3) 

set(gca,'FontSize',24) 

xlabel ('z') 

ylabel('sin/cos') 

xlim([0 1]) 

legend('sin(theta)','cos(theta) ' , 'brk' ) 
title('does \theta = \theta_o') 
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figure (2) 

set(gca,'FontSize',24) 

plot3(x,y,z) 

xlabel( ' x') 

ylabel( ' y') 

zlabel ('z') 

grid 


figure (3) 

subplot (2,2,1) 

plot(z,frac_ep) 

set(gca,'FontSize',24) 

xlabel (' z (m )') 

ylabel('frac ep') 

xlim([0 max(z)]) 

subplot (2,2,2) 

plot(z,frac_a) 

set(gca,'FontSize',24) 

xlabel (' z (m )') 

ylabel (' frac_a') 

xlim([0 max(z)]) 

subplot (2,2,3) 

plot(z,frac_b) 

set(gca,'FontSize',24) 

xlabel (' z (m )') 

ylabel (' frac_b') 

xlim([0 max(z)]) 

subplot (2,2,4) 

plot(z,frac_eta) 

set(gca,'FontSize',24) 
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xlabel(' z (m) ' ) 
ylabel(' frac eta') 
xlim([0 max(z)]) 

figure(4) 

plot(z,brk,z,brk2) 
legend('brk','brk2 ' ) 
xlim([0 max(z)]) 

figure(5) 

plot(z,stheta,z,stheta2) 
legend(' sthetal ',' stheta2') 
xlim([0 max(z)]) 

figure (6) 
plot(z,theta) 
xlim ( [0 0.1]) 
ylim ( [0 2*pi]) 

figure (7) 

plot(z,eta,z,etap,z,etam) 
legend (' eta ',' etap ',' etam') 
xlim([0 max(z)]) 

figure(8) 

plot (z,ggamma ,' LineWidth' , 3) 

set(gca,'FontSize' ,26) 

set (gcf, 'color'' , 'w' ) ; 

xlim([0 max(z)]) 

xlabel('z (m)') 

ylabel('\gamma' ) 

grid 
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figure (9) 

set (gcf, 'color'' , ' w' ) ; 
set(gca,'FontSize',32); 

plot(z_anal,ggamma*mc2_ev/lE9,'k','LineWidth' , 5) 
hold on 

plot(z,ggamma*mc2_ev/lE9,'—r','LineWidth',2) 
hold on 

%plot(z10,ggammal0*mc2_ev/IE9,'g','LineWidth',3) 

xlabel('z (m)') 

ylabel (' \gamma mc / '2 (GeV) ' ) 

legend('Analytic','Simulation') 

xlim([0 max(z)]) 

ylim([min(ggamma)*mc2_ev/1E9*0.7 1.l*max(ggamma)*mc2_ev/IE9]) 
grid 
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